HIDDEN SYMMETRIES AND DECAY FOR THE WAVE 
EQUATION ON THE KERR SPACETIME 

LARS ANDERSSONt AND PIETER BLUE* 

Abstract. Energy and decay estimates for the wave equation on the exterior 
region of slowly rotating Kerr spacetimes are proved. The method used is a 
generalisation of the vector-field method, which allows the use of higher-order 
symmetry operators. In particular, our method makes use of the second-order 
Carter operator, which is a hidden symmetry in the sense that it does not 
correspond to a Killing symmetry of the spacetime. 

The main result gives, in stationary regions, an almost inverse linear decay 
rate and the corresponding decay rate at the event horizon and null infinity. 
Except for the small loss in the decay rate, this generalizes the known decay 
results on the exterior region of the Schwarzschild spacetime. 
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1. Introduction 

In this paper we prove boundedness and decay for solutions of the covariant wave 
equation 

in the exterior region of the Kerr spacetime. In Boyer-Lindquist coordinates, the 
exterior region is given by (t, r, 0, S M x (r+, oo) x S*^ with the Lorentz metric 

+ '^#= + §d^' + Ed9^ (1,1) 

where r+ = M + \J Af ^ — a2 

A^r'^ -2Mr + a^, S = + cos^ 6*, Il = {r^ + a^f-a^sm^eA. 

For < |a| < M, the Kerr metric describes a rotating black hole, with mass M and 
angular momentum Ma, and with horizon located at r = r-|_. The Schwarzschild 
spacetime is the subcase with a = 0. The exterior region is globally hyperbolic, 
with the surfaces of constant t, St, as Cauchy surfaces. Thus, the wave equation is 
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well posed in the exterior region, even though the Kerr spacetime can be extended. 
We consider initial data on the hypersurface Eq • 

The Kerr black-hole spacetime is expected to be the unique, stationary, asymp- 
totically flat, vacuum spacetime containing a nondegenerate Killing horizon [2]. 
Further, motivated by considerations including the weak cosmic censorship conjec- 
ture, the Kerr black hole is expected to be the asymptotic limit of the evolution of 
asymptotically flat, vacuum data in general relativity. An important step towards 
establishing the validity of this scenario is to prove the stability of the Kerr solution, 
i.e. to show that vacuum spacetimes evolving from data which represent a small 
perturbation of Kerr initial data asymptotically approach a Kerr solution. Decay 
for the scalar wave equation on the Kerr background is an important model problem 
for stability, and due to the importance of this application, the wave equation on 
black-hole backgrounds has been actively studied in the last decade. 

An essential tool in the analysis of both linear and nonlinear Lagrangian field 
equations is the use of Noetherian currents associated to Killing or conformal sym- 
metries of the background spacetime. In the relativistic setting, we interpret these 
currents as momenta. A method, which may be referred to as the vector-field 
method, based on the systematic use of such currents, has been developed and 
has played an essential role in the proof of the nonlinear stability of Minkowski 
space [12], which built on earlier vector-field based estimates for the decay rates 
of solutions to linear and nonlinear wave equations p9] and to Maxwell's equa- 
tions and the spin-2 field equations [TT] . Generalizations of these ideas have played 
a central role in recent work concerning the wave and Maxwell equations on the 
Schwarzschild spacetime, and the wave equation on the Kerr spacetime. For these 
non-flat background spacetimes, the lack of symmetries presents an important new 
problem. 

The 10 dimensional group of isometries of the Minkowski space is broken to a 4 
dimensional group for the Schwarzschild spacetime, generated by dt and the spatial 
rotations. Further, the Schwarzschild spacetime contains orbiting null geodesies, 
located at the photon sphere, the hypersurface with r = 3M. As high frequency 
waves can track null geodesies for long times, an analysis of this feature is an 
essential step in a proof of decay for the wave equation. 

In the general Kerr spacetime, with a 7^ 0, which we consider in this paper, there 
are only two Killing flelds, dt and d^. In addition, in studying the wave equation on 
the Kerr spacetime, one encounters several new phenomena which are not present 
in the Schwarzschild case. There is an ergo-region outside the horizon, where the 
stationary Killing field dt fails to be timelike. Thus, the Kerr spacetime admits no 
positive definite, conserved energy for the wave equation. Further, the orbiting null 
geodesies in Kerr fill an open region in spacetime. The lack of symmetries of the 
Kerr spacetime is compensated for by the presence of a fundamentally new feature, 
which we make essential use of in this paper, a hidden symmetry. 

By a hidden symmetry we mean an operator which commutes with the wave op- 
erator, not associated to a Killing vector field, but rather to a second-rank Killing 
tensor. For the Kerr spacetime, the Killing tensor and the related conserved quan- 
tity found by Carter ^ provides, via the associated second-order Carter operator, 
a hidden symmetry. The existence of the two Killing vectors and the Killing ten- 
sor imply the separability of many important equations on the Kerr spacetime, 
including the wave equation. 

One of the fundamentally new ideas introduced in this work is a generalization 
of the vector-field method which allows the use of not only Killing symmetries but 
also hidden symmetries in the construction of suitable Noetherian currents for the 
analysis of Lagrangian field equations. This allows us, in contrast to other recent 
work on the wave equation on Kerr, to carry out our proof of uniform boundedness 
and decay results exclusively in physical space, using only the coordinate functions 
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and differential operators. This technique almost eliminates the need for methods 
involving separation of variables or Fourier analysis |^ 

To state our main results, we use the tortoise coordinate r*, defined by 

dr r^~2Mr + a^ 

r,(3Af) = 0, 



dr^ + 
and the almost null coordinates u± given by 

M+ = i + n, u^=t — r^. 

Through out this paper, we will take M fixed, so that, in a statement about the 
existence of some a for which an estimate holds for |a| < a, it is understood that 
the upper bound a depends on M . Our main results are: 

Theorem 1.1 (Uniformly bounded, positive energy). There are positive constants 
Ci and di, and a nonnegative quadratic form on each hypersurface of constant t, 
Et^ [4']{t), such that, if \a\ < di and ^ : M x (r+, oo) x 5^ ^ M is a solution of the 
wave equation, V"Vq^ = 0, then Vt 

Theorem 1.2 (Decay estimates). There are positive constants C2, C2, and d2, and 
there is a nonnegative quadratic form on each hypersurface of constant t, ||-0|p(t), 
such that, if\a\ < 0,2 andip ■ Rx (r_|_, 00) x.?^ ^ R is a solution of the wave equation, 
V^Vq^ = 0, then there are the following decay estimates \ft > 0, {9, (j)) £ : 

(1) Decay in stationary regions: Vr G (3Af, 4M); 

m,r,e,cl>)\ <C2max{l,i}"'+'^^l°l||V'||(0). 

(2) Near decay: Vr < MI: 

m,r,e,cl>)\ <C2max{l,u+}-i+^^l'^l||V;||(0). 

(3) Far decay: Vr with r > AM and r < t: 

|V'(t,r,g,^)| <C2r-^max{l,^_}g^H ^+ ^ ' |i^||(o). 

\u+ maxjl, U-\ J 

In particular, for t/2 < r < t: 

r, e, cj))\ <C2r-' max{l, |1^||(0). 

(4) Decay near spatial infinity: Vr with r > t: 

mt,r,e,cl>)\ <C72r-imax{l,-ii_}-i/2|lV||(0). 

Theorem 1 1.1 1 is the conclusion of section |3] and is given in theorem |3.14| Theorem 

Except for the loss 
are the same as the 



1.2 follows from the conclusions of theorems 15.1 5.2 and 



5.4 



1.2 



in the exponent of C2|a|, the estimates stated in theorem^ 

results proven using vector-field techniques in the Schwarzschild spacetimej^ The 
decay along null infinity is the same (modulo the loss in the exponent) as can be 
obtained in Minkowski space from initial data on t = which decays like r~^/^. 
This is roughly the decay rate we require for the initial data in Kerr. Note that the 
Kerr spacetime has a discrete, time reversal symmetry (i, 0) i-^ (— — 0), so that 
the results of theorems 11.11 and 1 1 . 21 can also be reversed to i < 0. 



^We do not use separability or the Fourier transform (which are essentially equivalent to each 
other) in the t coordinate. In the proof of lemma |3.13[ we need to use separability in the 
coordinate to control the axially symmetric component of the solution, but are otherwise able to 
avoid using separability (or the Fourier transform) in </>. We never use separability in the r and 
8 coordinates, since this is only possible once separation in the t and </> coordinates has already 
been performed. 

■^Recently, better decay estimates have been proven in the Schwarzschild case, giving, for any 
positive S, decay rates of J— 3/2+5 ^^^^ tt^^'^^^* in stationary and near regions respectively I33| . 
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The norm Ht/iIKO) in theorem 1.2 is bounded if ■0 is smooth on the hypersurface 
So and if ^p and its first nine derivatives (with respect to the Boyer-Lindquist 
coordinates) decay like 7--3/2+5 ^~5/2+s respectively for some positive S as r —>■ 
oo. By "smooth" we mean C°° with respect to local coordinates. As r ^ this 
is not the same as being smooth with respect to the Boyer-Lindquist coordinates, 
since they degenerate here. The quadratic form Et^ [V'KO) in theorem ] l.l| is bounded 



when IIV'IP(O) in theorem 1.2 is bounded. 

We briefly comment on other related work. Estimates for the decay rate of so- 
lutions to the wave equation have been proven in the subcase of the Schwarzschild 
spacetime, where a = 0. Birkhoff's theorem states that the Schwarzschild space- 
time is the unique spherically symmetric, vacuum spacetime solution of Einstein's 
equation. For the coupled Einstein and scalar wave system, a decay rate and non- 
linear stability of the Schwarzschild solution have been proven in the spherically 
symmetric setting [IB] . 

As mentioned earlier, for the wave equation without a symmetry assumption 
but on a fixed background spacetime, the case of the linear wave equation on the 
Schwarzschild spacetime is significantly simpler than the corresponding case in the 
Kerr spacetime, since the dt Killing vector is timelike in the entire exterior region 
and generates a conserved positive energy, there is the full SO{3) group of rotation 
symmetries available to generate higher energies, and the orbiting null geodesies 
are restricted to r = 3M. The first two of these properties imply that solutions 
remain bounded. Following the introduction of a Morawetz vector field and of the 
equivalent of an almost conformal vector field to the Schwarzschild spacetime |31j , 
decay estimates for the wave equation were proven , proven with a weaker decay 
rate but less regularity loss , and proven separately with a stronger estimate near 
the event horizon [18]. These were extended to Strichartz estimates for the wave 
equation [34] and decay estimates for Maxwell's equation [4] . The Morawetz vector 
field which made these estimates possible was centred about the orbiting geodesies 
aX r — 3M. This construction of a classical vector field fails when a 7^ 0, since there 
are orbiting geodesies filling an open set in spacetime. 

Recently, new Morawetz vector fields with coefficients that depend on both space- 
time position and on Fourier operators have been introduced to construct a uni- 
formly bounded, positive energy and to show that solutions to the Kerr wave equa- 
tion remain uniformly bounded [T7] ST] . These might reasonably be called Fourier- 
analytic, pseudodifferential, microlocal, or phase-space techniques, since the Fourier 
operators represent coordinates in momentum space in contrast to spacetime co- 
ordinates in physical or configuration space. These results include a form of weak 
decay, since the Morawetz estimate implies integrability of the local energy. In |16j , 
it has been announced that these Fourier-analytic techniques can be extended to 
decay results very similar to ours in theorem |1. 2 [ Like our work, these require that 
\a\ is very small relative to M . We compare this work with our own in more detail 
in section [oj Very recently, a decay rate estimate of for all \a\ < M has been 
announced |15|. and a proof of decay has been given [4^. Both of these take a 
local energy estimate, like the Morawetz estimate in this paper, as an assumption. 
The techniques in this paper remain the only ones to provide such an estimate with- 
out Fourier transforms in time. Fourier- analytic vector fields were used previously 
to prove Mourre estimates, which are similar to Morawetz estimates, in the proof 
of scattering for the Klein-Gordon equation ^26) and the Dirac equation [27! ■ 

Finally, we recall that decay for the wave equation has previously been obtained 
[22 ] from an explicit representation of solutions ^l] using the complete separability 
of the Kerr wave equation. These decay results are of the form limt_>oo \^p{t, r, 0, (f))\ 
0, where ip{t,r,9,(p) — i/ji^ (t, r, ^)e*^-'^ or where tp is made up of a finite number 
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of azimuthal modes of this form. Decay rates have been obtained from this sepa- 
rabihty method for solutions to the Dirac equation |20j and spherically symmetric 
solutions to the wave equation when a = [SU] ■ The decay without rate results for 
the wave equation built on the earlier result that there axe no exponentially growing 
modes |45] and applies for all \a\ e [0,M]. 

1.1. Hidden symmetries and the vector-field method. In the classical vector- 
field method, one seeks to control the solution of field equations by making use of 
energy fluxes and deformation terms calculated for suitably chosen vector fields. The 
vector fields used are often (approximate) conformal symmetries of the background 
spacetime. Higher-order estimates are achieved by Lie differentiating the field along 
further vector fields. An important advantage of the vector-field method is that one 
works entirely in terms of quantities in physical space. 

A particularly clear example of the vector-field method was the use of the Lorentz 
group in the Minkowski spacetime, M^"*"^, to construct norms which could be used 
in the Klainerman-Sobolev inequality and to use this to prove the well-posedness of 
nonlinear wave equations [29 . The term "vector-field method" seems to have been 
introduced relatively recently, especially to describe generalisations of the work in 
to situations where one lacks the full Lorentz group of symmetries. Applying 
this terminology retroactively, we would now describe the early uses of the radial 
Morawetz vector-field [36] and the conformal vector-field [25] as applications of the 
vector-field method. This may have previously been refered to as the method of 
multipliers or the Euler-Lagrange method, although these terms can also be applied 
to more general techniques. 

A central result in mathematical relativity, and perhaps the most important ap- 
plication of the vector-field method, was the proof of the nonlinear stability of the 
Minkowski spacetime [T2] . There had also been earlier work on the stability of the 
Minkowski spacetime, but this required hyperboloidal initial data |,23 . The mon- 
umental proof of nonlinear stability built upon previous vector-field estimates for 
linear and nonlinear wave equations |29J, and for the Maxwell and spin-2 field equa- 
tions [11 , which are better models for Einstein's equations. This partly motivates 
our work on the linear wave equation in the Kerr spacetime using generalisations 
of the vector-field method. Since the original proof of nonlinear stability for the 
Minkowski spacetime, a simpler proof has been developed, but this also makes use 
of the vector- field method [321. 




As mentioned above, in the Kerr spacetime, the lack of symmetries, as well as 
the complicated nature of the orbiting null geodesies, makes it impossible to derive 
the required estimates using only classical vector fields. In this section, we outline 
a generalization of the vector-field method which allows us to take advantage of the 
presence of hidden symmetries in the Kerr spacetime. In particular, we consider 
energies based on operators of order greater than one, rather than just vector fields. 

Let Og — V"Vq. In the discussion here, we consider the scalar wave equation 
Ogip = 0, but much of the discussion applies equally to general field equations 
derived from a quadratic action. We define a symmetry operator to be a differential 
operator S such that if Ogtjj = 0, then also OgStp = 0. The set of symmetry 
operators is closed under scalar multiplication, addition, and composition, and each 
symmetry operator has a well-defined order as a differential operator. Thus, the set 
of symmetry operators forms a graded algebra. Given a set of generators of the set 
of symmetries, we can consider the subset consisting of generators of order n. We 
denote this subset of the generators of the symmetry operators by S„ and denote 
the elements of S„ with an underlined index, e.g. Sa G §«. 

If X is a conformal Killing field, then the operator £x generated by Lie dif- 
ferentiation with respect to X is clearly a symmetry operator. We take a hidden 
symmetry to be a symmetry operator which is not in the algebra generated by 
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the Killing vector fields. Since the Minkowski spacetime saturates the Delong- 
Takeuchi-Thompson inequality, there are no hidden symmetries [lOj . Similarly, in 
the Schwarzschild spacetime, there are no hidden symmetries [5]. In contrast, in 
the Kerr spacetime, it is well-known that there is Carter's Killing 2 tensor, and that 
this generates a hidden symmetry [8l HI] . 

The energy-momentum tensor for the wave equation is 

T[ij]al3 = V„^V/3V^ - ^gal3 (V^VV'^V') • (1-2) 

The momentum associated with a vector field X and the energy associated with a 
vector field X and evaluated on a hypersurface E are 



S 

where &q is integration with respect to the surface volume induced by g on S. In 
the following, unless there is room for confusion, we will drop reference to ip in the 
notation for momentum and energy. When the spacetime is foliated by surfaces of 
constant time, we will denote these surfaces by St and typically denote the energy 
on such a surface by -Ex(i) = -Ex(S(). 



The energy momentum tensor (1.2 1 satisfies the dominant energy condition, and 
hence for X timelike, the energy induced on a hypersurface with a timelike normal 
(i.e. a spacelike hypersurface) is positive definite. The energy conservation law 
takes the form 

i?x(S2)-£;x(Si)= / (V„P5)y^d4x, 
Jn 

where is the region enclosed between Ei and S2. This is often referred to as the 
deformation formula. Energy estimates are often performed by controlling the bulk 
(also called deformation) terms V^Px. However, for the Morawetz estimate (e.g. 
inequality ( |3.6| )), one makes use of the sign of the bulk term itself. 

By estimating higher-order energies one may, via Sobolev estimates, get pointwise 
control of the fields. Higher-order energies may be defined by using symmetries. If 
for < i < n, there is a collection of order-n differential operators. Si, then we can 
define the higher-order energy (of order n) for a vector field X to be 

n 

i?x,„+i[V'](S) = 5]^i?x[5V^](S). 

Since the energy-momentum tensor is quadratic in if), we can define a bilinear 
form of the energy momentum by 

It is convenient to define an index version of the bilinear energy momentum, with 
respect to a set of symmetry operators {Sg} by 

Given a double-indexed collection of vector fields, {X—}, we define the associated 
generalized momentum and energy to be 



i?x-M(S) = / Px-Mad?7" 



In practice it is convenient to consider momenta with lower-order terms, designed 
to improve certain deformation terms in V^Px- For a scalar function, q { (34j . but 
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previously appearing in |14j). or a double-indexed collection of functions, q— , the 
associated momenta are defined to be 

P,[V']a='Z(V„V')^'-^(aag)V', 

For a pair consisting of a vector field and a scalar function, (X,^/), the associated 
momentum is defined to be the sum of the momenta associated with the vector field 
and the scalar. For a pair of collections, (X^^,g^^), again the momentum is defined 
to be the sum of the momenta. In all cases, the energy on a hypersurface is given by 
the fiux, defined with respect to the momentum vector, through the hypersurface. 

It is important to point out, as we show in lemma [2TT| that the deformation terms 
for the generalized momenta are computationally not much more difficult to handle 
than the classical ones. As for the classical momenta and energies, in defining the 
generalised vector fields, momenta, and energies as outlined above, one is interested 
in getting positive definiteness of the energies or bulk terms. Here, an additional 



subtlety arises. Namely, in the Morawetz estimate presented in equation (3.6 1, 
one achieves positive definiteness only modulo boundary terms. We generate these 
boundary terms when we integrate by parts to use the formal self-adjointness of the 
second-order symmetry operators. These boundary terms can then be controlled 
by the energy. The presence of these boundary terms is a completely new feature 
compared to the classical energies and momenta. 

1.2. Symmetries and null geodesies of Kerr. For any geodesic, the quantity 
5q/37"7^ is a constant of the motion. Given a Killing field ^, the quantity — 
gapi"^^ is an additional conserved quantity. In Kerr, we have the Killing fields dt 
and dif, with the associated constants of motion pt and p^. For a timelike or null 
geodesic, these correspond to the energy and the angular momentum of a particle 
or photon with world line 7 and are denoted E and L^. 

More generally, if the spacetime admits a Killing fc-tensor, i.e. a symmetric 
tensor Kai - ok which solves the Killing equation (aiKa^-'-ak+i) = 0; then K = 
^ai---afc7"^ ■ • • 7"'° is a conserved quantity. In the particular case of Killing 2- 
tensors, which is the only case we are interested in here, there is associated to the 
Killing tensor a symmetry operator K — W aK'^'^W p, such that [iC, Dg] — [51 HI]. 
Since the commutator is zero, this operator is clearly a symmetry in the slightly 
weaker sense defined in the previous section. 

In the Kerr spacetime. Carter's Killing 2-tensor, provides a fourth constant of 
the motion Q = Qa/si^'j^ ■ For a null geodesic, we have 

O COS i-) O . O y\ o 

^ _ „^ , -'^ ' a sm^ '^"■^ 



sm 

A similar expression exists for timelike or spacelike geodesies. Any linear combi- 
nation of E"^ , ELz, and can be added to Q to give an alternate choice for the 
fourth constant of the motion. The form we have chosen is uncommon, but useful 
for our purposes because it is nonnegative. 

As was demonstrated by Carter, the presence of the extra conserved quantity 
allows one to separate the equations of geodesic motion. Of most interest to us is 
the equation for the r coordinate of null geodesies, 

S'f^) =-n{r-M,a-E,Lz,Q), (1.3) 



where 



n{r;M,a;E,Lz,Q) = 



- (r^ + a^)^E^ - iaMrELz + (A - a^)Ll + AQ. (1.4) 
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One finds that orbiting null geodesies, i.e. ones which do not fall into the black hole 
or escape to infinity, must have orbits with constant r. The r values allowing orbiting 
null geodesies are solutions to the equations 7?. = 0, &lZ/dr = 0. The solutions to 
this system in the exterior region turn out to be unstable, which corresponds with 
our conventions to d^TZ/dr'^ < 0. 

In the Schwarzschild case, i.e. for a — 0, there are only orbits on the sphere 
at r = 3Af , which is called the photon sphere. For nonzero a, the orbiting null 
geodesies fill up an open region in spacetime which we shall also refer to as the 
photon sphere in the Kerr case. As a — *■ 0, the photon sphere tends to r = 3M. 
There are many descriptions of the Kerr spacetime and its geodesies, including 

In Boyer-Lindquist coordinates, the d'Alembertian Dg — V"Vq, takes the form 
ag = ^ (drAdr + ^n{r;M, a; dt,d^,Q)) , (1.5) 



where TZ is given by (1.4 1 with the conserved quantities E,Lz,Q replaced by their 
corresponding operators dt,d^, and the second-order Carter operator]^ Qi 



1 . cos2 6' 

^ty^-ru, Dill uu^ 



ede sin Ode + '^^^^^dl + sin^ Qd}, 



sm sm 

7^(r; M, a; 9*, 9^, Q) = - (r^ + a^fd} - AaMrdtd^ + (A - a^)dl + AQ. (1.6) 

We have used some unusual sign conventions in defining TZ to avoid factors of i 
when replacing the constants of motion by differential operators. 

It is clear from the above that dt, d^, and Q are symmetry operators for the 
wave equation on Kerr. We denote the set of order-n generators of the symmetry 
algebra generated by these operators by 

= {drd2*Q''Q\nt + n^ + 2nQ = n;nt,n^,nQ e N}. (1.7) 

In particular, 

§o={Id}, §i={dt,d^}. 

Of particular importance in our analysis will be the set of second-order symmetry 
operators, 

S2 = {dldtd^,dlQ} = {Sa}, 

and underlined indices always refer to the index in this set. 

The function TZ is polynomial in its last three variables, so TZ{r; M, a; dt,d^, Q) is 
well defined. Furthermore, it can be written as a linear combination of the second- 
order symmetries with coefficients which are rational in r, M, and a, 

TZ{r-M,a;dt,d^,Q) ^TZ^Sa- 



1.3. Strategy of proof and further results. Recall from earlier in the introduc- 
tion that there are three major problems in the Kerr spacetime: 

(1) No positive, conserved energy: There is no timelike. Killing vector. In 
particular, the vector field dt, which is Killing, is only timelike outside the 
ergosphere, r > M + ^/AP — cos-^ 6. 

(2) Lack of sufficient classical symmetries: The higher energies generated the 
Lie derivatives in the dt and 9^ directions do not control enough directions 
to control Sobolev norms or the function. 



Since the Carter operator Q and the Carter constant are closely related, we use the same 
notation for both. 
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(3) Complicated trapping: There are null geodesies which orbit the black hole, 
in the sense that they neither escape to null infinity nor enter the black 
hole. Since solutions to the wave equation can follow null geodesies for 
an arbitrarily long time, this presents an obstacle to decay. Furthermore, 
there are orbiting geodesies occuring over a range of r in the Kerr spacetime 
(with \a\ > 0), which makes the situation more complicated than in the 
Schwarzschild spacetime (a = 0), where there are only orbiting geodesies at 
r = 3M. 

To overcome the first problem, we first observe that the vector dt is timelike for 
sufficiently large r; that, if 

a 



rl 



denotes the angular velocity of the horizon, then the vector dt + iond^ is null on 
the horizon and timelike for sufficiently small r > r+ ; that the regions where dt and 
dt + i^Hdr/i are timelike overlap when \a\ is sufficiently small; and that both dt and 
dt + i^ndcf, are Killing. Thus, if we let 

T^^dt + x^Hd^, (1.8) 

where x is identically 1 for r < for some constant r^, identically for r > 
+ M, and smoothly decreases on [r^, + M], then, for sufficiently small a, this 
vector-field will be timelike everywhere and will be Killing outside the fixed region 
r G ['"x'^x + Thus, to prove the boundedness of this positive energy, it will 
be sufficient to control the behaviour of solutions in this fixed region through a 
Morawetz estimate. 

To overcome problem ([2| , we note that the second-order operator Q is a symme- 
try and is a weakly elliptic operator. Using Q, 9^, and df as symmetries to generate 
higher energies, we can control energies of the the spherical Laplacian of tp. These 
control Sobolev norms which are sufficiently strong to control |V-'P- 

To handle the complicated trapping, we will use our extension of the vector-field 
method to include hidden symmetries. To construct a Morawetz multiplier, we 
would like to construct a vector field with a weight that changes sign at the orbiting 
geodesic, but this is not possible using a classical vector-field. If we introduce 
C — C-Sa = df + d^ + Q to give us an elliptic operator and an extra, free, underlined 
index, we can take as our collection of Morawetz vector fields 



= - zwn'^'-^CP^ dr, 



with z and w smooth, positive functions to be chosen. Applying the analogy of the 
deformation formula, the difference between the energies on one hypersurface and 
another is 



Ignoring several distracting details, the deformation is of the form 

]^zn"i-n"^c''^{do.SamdpSbi;) 



4 

X. 



-(dr^drZ{drW'R"^))C\Sa^){Sb'ip). 
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In the first line, one factor of TV arises from the wave equation, and the other from 
our choice of the Morawetz vector field A—, which allows us to construct a perfect 
square to obtain positivity. In the second line, the term involves two derivatives of 
—iZ. Near the photon orbits, the convexity properties of 7?., which ensured that the 
orbits are unstable, ensure that this term is positive. We are free to choose z and w 
to get positivity away from the photon orbits. The third term is lower order, since 
it involves fewer derivatives. 

For small |a|, with v denoting terms of the form Sai^, and with our choices of z 
and w, the sum of the second and third terms is of the form 

M [-.a^^frvf + v^) (1.9) 

with small perturbations on the coefficients. The coefficient on is positive outside 
a compact interval in (r+,oo). As shown in [5], it is sufficient to prove a Hardy 
estimate which bounds the quadratic form in ( 1.9 1 from below by a sum of positive 
weights times {drv)^ and v'^. 

The positive terms arising from the deformation of A— dominate the deformation 
terms (with extra derivatives) terms arising in the failure of to be Killing. 
(In fact, the terms in the Morawetz estimate only control the second and third 
derivatives of ip, where as the deformation terms from the third-order energy 
also involve the first derivatives and undifferentiated factors of ip. It is at this 
point, in the proof of lemma 3.13 where we are forced to make a decomposition 
in harmonics oi dip, i.e. to separate variables in 0, to obtain additional control on 
the rotationally symmetric components.) On the other hand, the energy associated 
with A— is dominated by the (third-order) energy associated with T^. Since there 
is a factor of \a\ on the T^^, deformation terms, we have a small parameter, which 
allows us to close the boot-strap argument in which the energy is controlled 
by the integral of its deformation, which is controlled by the integral of the A— 
deformation, which is controlled by the A— energy, which is finally controlled by 
the energy. This allows us to establish theorem |1.1[ 

To prove decay, we introduce the vector field 

K = (t^ + rl + 1)9, + 2tr., + 5... 

In the Minkowski spacetime, the corresponding vector field is a conformal Killing 
vector field and generates a positive, conserved energy, sometimes called the confor- 
mal energy or conformal charge. (Confusingly, both A— and K are called Morawetz 
vector fields. We will only use the term Morawetz vector field to refer to A—.) In the 
Schwarzschild spacetime, it is now well known that the corresponding vector field 
has a deformation tensor which is favourable outside a compact region in (r+,cx)), 
and that growth in the K energy can only occur as a result of the wave remaining 
inside the compact interval for long periods of time. Since the Morawetz estimate 
rules this out, after another boot-strap argument, it is possible to bound the K 
energy. After changing variables, we are able to mimic most of the argument from 
the Schwarzschild spacetime. We need to introduce some lower-order terms, (/k, to 
remove the worst terms in the deformation expression, but are still left with terms, 
in the K deformation, of the form 

at^{dri^){d,j,il)). 

since these involve a quadratic expression in the derivatives times a factor of t^, 
these are of the same order as the K energy itself. Even after applying the Morawetz 
estimate, one is left with an estimate of the form 

E-K{t2) — Ei^{ti) < C\a\ / dt + (more easily controlled terms), 
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where Eii{t) denotes the energy evaluated on the hypersurface {t} x (r+,oo) x S^. 
Although one cannot obtain a uniform bound on the K energy from this, if the 
surfaces are of constant t, then the growth cannot be faster than Since, in 

regions of fixed r, the K energy is like the T^^ energy times t^, this means that the 
energy density in regions of fixed r (away from the horizon) decays like i-2+c|a|^ 
Since the energy density is quadratic, if we use the third-order K energy, we can 
control the solution tp in regions of fixed r by ^^i+ci"!. which proves the first part 
of theorem 11.21 

We also prove decay near to and far from the black hole. In both cases, we drop 
a null geodesic (or almost-null curves) from the point q at {t,r^,,9,<j)) to a point 
p dX r = 3M, where we already have decay by theorem |1.2[ and then estimate 
\^{q) — "0(^)1, where ■0 = (r^ -t- a^Y^'^ijj. The point p at r = 3M has a t coordinate 
given by the m+ or m_ coordinate of q in the near or far cases respectively. In 
the near case, following [H [M], we consider the region in the (u+,w_) plane and 
enclosed by p, q, curves of constant u+ and of constant m_, and t = 0. We apply 
Stokes' theorem in this region to the one-form {d-tfj)du- and estimate the bulk 
term using the Morawetz and K estimates. This allows us to control {tpiq) — ip{p)\ 
as one piece of the boundary term. In the far case, we drop an almost-null curve 
from q to the surface r = t/2, and we call this intersection point p. Although p is 
not in a stationary region, we are still able to show sufficiently strong decay for ip{p) 
directly from the K energy bound. We evolve the solution from the hypersurface 
t = to the hypersurface of constant t through p. Since the deformation of the K 
energy is easily controlled when r > t/2, we can deform the hypersurface so that 
it becomes almost null and passes through q. We then integrate the derivative of 
ip along an almost null curve in this hypersurface to estimate {ipiq) — '<P{p)\- By 
applying the Cauchy-Schwarz inequality, we are able to control the integral of this 
derivative by the product of the K energy on the hypersurface and by the desired 
decaying factor. The contribution from the endpoint p at r = 3M decays much 



faster. These estimates give the remaining parts of theorem 1.2 



The small \a\ condition which we impose is significantly stronger than the con- 
dition that \a\ < M which implies the existence of a black hole and which might 
be ideally imposed. There are several fundamental and technical reasons for this 
small a condition. Perhaps most importantly, the construction of relies on there 
being a region where both dt + X'^ndip and dt are time like in which to perform the 
blending. When \a\ is sufficiently large, but still smaller than M, there is no such 
overlapping region, so this particular construction fails. In addition, we use the 
assumption on the smallness of |a| to close the bounded energy argument and 
the K energy growth argument. If \a\ is not small relative to the absolute constants 
appearing in those estimates, it would not be possible to close the boot strap. A 
clear technical obstacle is that, in the proof of the Morawetz estimate, we perturb 
the Hardy estimate in ( 1.9 1. If \a\ were too large, the perturbation argument would 
fail, and our numerical investigation suggests that when \a\ is larger than about 
.9M, there are no longer positive solutions of the associated ODE, which we use to 
prove the estimate. These obstacles are the most fundamental obstacles to extend- 
ing the range of |a|, but there are also numerous other, technical estimates in which 
we have made use of the smallness of \a\. 

Having summarized our method, we will now compare it with methods used in 
recent, related work. Recently, others have constructed a bounded energy |17l 141] . 
To make a comparison, we point to several features which they share but which are 
different from those in our approach. 

Since our energy is based on T^, which becomes null on the event horizon, the 
energy we control has a weight which vanishes linearly at r — r^. The other works 
make use of the horizon-penetrating vector field, first introduced in |18j . This is 
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denoted Y [17 or X2 |41j . By combining the horizon-penetrating vector field with 
the equivalent of T^, they are able to construct a timelike vector field and an energy 
which do not degenerate near the event horizon. This is clearly advantageous. Since 
the properties of the Y or X2 vecter field are quite stable [16] , it can be used in a 
separate step, which we omit, following our result. 

Neither [T7| nor [IT] use Q to generate higher energies. Away from the event 
horizon, they use the symmetries 9f , 9t9^, and 9^ and the fact that ip satisfies the 
wave equation. Near the event horizon, they generate higher energies using dt and 
a horizon-penetrating, radial vector field (e.g. Y in [T7]). This is possible because 
of a favourable sign in the error terms arising from the failure of the radial vector 
field to be a symmetry. 

Both of [TT] m] perform a Fourier transform in the t and (f> variables to con- 
struct a pseudodifferential Morawetz multiplier, which we have avoided in favour of 
differential operators. 

Less importantly, both avoid surfaces of constant t in favour of surfaces and 
coordinates which go through the event horizon. Since vector-field arguments can 
be deformed from one surface to another, this is a minor difference, although, the 
presence of the functions, q, slightly complicates this. Although all known Morawetz 
arguments have, in some sense, a troublesome lower-order term, |171 141] use a 
different construction so that they can use positivity arising from Y or X2, instead 



of the Hardy estimate we use to control the negativity in (1.9). 

In [16], a decay rate of t^i+*(°) is stated with 5{a) ^ as a ^ and with the 
corresponding decay rates along the event horizon and null infinity. The detailed 
outline of a proof focuses on using a vector field K, which is very similar to the 
one used previously in the Schwarzschild spacetime and, hence, to the one we use 
in this paper. 

The structure of this paper is as follows: In section [2] we provide some further 
notation which we use in this paper. In section [3j the main argument of this paper, 
we expand the energy associated with T,^ and prove the Morawetz estimate using 
the symmetry-indexed vector fields. This is followed by the K argument to prove 
local energy decay in section |4j and then by the decay estimate for 'ip itself in section 

m 

2. Notation and preliminaries 

In this section, we present some more notation and basic estimates which we will 
use through out the paper. 

To begin, we note that, in estimates, C is used to denote an absolute constant 
or a constant which depends only on M. The notation x < y means x < Cy, and 
the notation a; ~ y means x < y and y ^ x. Further, it is convenient to introduce 
the following notation 

d^/i ^fidddcj), fi ~ sin 9, 

d^ji —d^fidr, d^/i* ^d^/idr^, 

d^'/i ~d^ fidrdt, d"^//* ^d^ iidr^dt. 

2.1. Canonical analysis. The volume element for the Kerr metric in Boyer-Lindquist 
coordinates is 



v/^=Ssin0. 



It is convenient to consider instead of the covariant d'Alembertian Dg, the trans- 



formed d'Alembertian □ = SDg. Recalling (1.5 1, we can write □ in the form 



□ = drAdr + ^7^(r; M, a; dt,d^,Q). 
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This is the form of the d'Alembertian that we shall consider throughout this paper. 
Let 

g'^P = Y.gO'P ^ /i = sin6'. 
Then = is the Euler-Lagrange equation for the action 

s = J /:dV, 

with the Lagrangian 

The canonical energy-momentum tensor for S is 

and the momentum vector for a vector field X is given by 
or explicitely 

We have the following relations with the quantities introduced in section 

VaPf-x.,)^ = - (^^(x,,)[^]) , 

When it is clear from context what the arguments are, we will frequently write 
i?x W\ I -^x {t) , or £'x ■ Similarly, we will typically write 'P(x,q) for 'P(x,g) ["0] • Towards 
the end of this paper, we also need to evaluate energies on other hypersurfaces. In 
this case, we generate a three form by contracting the vector field 'P(x,(/) against 
the volume form generated by the coordinate one forms, and then integrate that 
3-form over the corresponding three-dimensional hypersurface. This is explained 
further in section 15.31 

The advantage of the canonical formalism just introduced is that without com- 
puting covariant derivatives, one can calculate the divergence of the momentum. 

Lemma 2.1. //X is a vector field, q is a function, and ifV(x.,q) the associated 
momentum relative to the Lagrangian and weight 

c = ]^{g"P{d^u){dpu) + vij^), 
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then the divergence of the corresponding momentum is 



+ {lCd-,4' + qi^) (^^5^/2^"%-^^ • (2.1) 

If 4^ is o, solution of [p.^^dafLG'^^dp — — then the last term in this formula 

vanishes. 

Similarly, ifX.— and q— are symmetric collections of double indexed vectors and 
scalars respectively and ^ is a solution of [fi~^dafJ.G'^^dp — ip = 0, then 



4 5a (/i^'fx^, 



.— ,9—, 



+ l^ido,fig''^di3q^){Sai^){SbjP). (2.2) 
2.2. The 3+1 decomposition. The normal to the surfaces of constant t is 

1/2 



/ n \ 



with uj± defined below. The lapse function is defined solely in terms of the foliation 
and its normal and is 

iV=(ng^9„t)"i 



Unfortunately, the t co-ordinates do not extend beyond the exterior of the black 
hole, so the vector field n^^ does not extend continuously beyond the exterior. 
We introduce the vector field 

T_L =dt + iL>±d^, 
2aMr 

This vector field has the properties that 

that Tj_ is timelike in the exterior, and that it extends continuously to the event 
horizon and the bifurcation sphere. In fact, it extends smoothly through the event 
horizon and the bifurcation sphere]^ This vector field extends to the null tangent 
vector on the event horizon and to axial rotation (with coefficient ujh) on the bi- 
furcation sphere. 



*The vector fields dt and are known to extend smoothly through the bifurcation sphere [27 
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2.3. Pointwise Norms. First, we introduce some notation for angular derivatives. 
We use y to denote angular partial derivatives and ^ for the spherical Laplacian 
in (9, 0) coordinates. For two vectors on the sphere, there is an inner product 
defined using this standard metric on the unit sphere. Extending this notation to 
differential operators in the standard way, for a function /, ^/ = (fi^f). 

We use Qi for the rotation vector-fields about the coordinate axes. With the 
exception of O3 ~ d^, these are not symmetries in Kerr. We use Oi — {Qi} to 
denote the set of these rotations, and we use Ti for {dt,Qi}. 

Given a set of differential operators, X, we use the notation 

\^\l=\x^\'=J2 

If no set is specified, simply an index, we mean 

n 
1=0 



where §„ is the set of generators of the order-n symmetries given in equation ( 1.7 1. 
We will refer to \ip\n as the n pointwise norm of ip. When the n is clear from context, 
we will simply refer to this as the norm of ip. 

Lemma 2.2 (Spherical Sobolev estimate using symmetries). There is a constant, 
C, such that V(0, G S"^ and all {t,r) G M x (r^^oo), if Tp is sufficiently smooth 
that the quantity on the right is bounded, then 



sup IV-I" 

(t,r)xS2 



'2 < c / i^i^dV 

J{t,r)xS^ 

If u is sufficiently smooth that the integral on the right is bounded. 



Proof. Recall that we use ^ to denote the spherical Laplacian 



-A- 



\m 



< 



-dgfide 
-de fide 



cot^ 0dl 
cot^ 0dl 



The absolute value of the spherical Laplacian of u can be estimated by 
By a standard, spherical, Sobolev estimate. 



(S2 



< 



Since the integrand on the right is bounded by 
a uniform constant in (t,r). 



I2, the desired estimate holds with 

□ 



In subsection |3.4| we also require the following operator and the associated weaker 
norms. 



Definition 2.3. Let 



C =d1 - 



16 



L. ANDERSSON AND P. BLUE 



and 

mL =e|5tVP + (1 + e) + ^l^tS^V'Pj + \m 

I^IL =^'\d'U'\' + (26 + e')\d^n\' + (1 + 2e)\dtM' + 17^'- 

We also introduce the homogeneous norms, generated from the previous norm by 
taking e= 1, 

IV'U.i- 

In theese norms, there are coefficients with lower-order terms in e, such as the 
(1 + e) in IV'b.e- The purpose of these lower-order terms is to permit exact equality 

in some estimates below. 

Lemma 2.4 (The C^C estimate). There is a positive constant C such that, for 
positive values of the parameter e, if \a\ < Ce and if ip smooth, then 

(A^)(£V) >mh (2.3) 

- -f- (/x(atV)(y9tV)) 

A* 

+ - {dt{i^2a^sm^e{dti^){dli;)) - d^{n2a^ sin'' 9 {dt^){d^dt^lj))) . 

Thus, (CciJj){Cip) is equal to the 2,e norm plus divergence terms. 

Proof. This follows by direct computation, but it is important to integrate by parts 
in t first. 

(£,V)(-C^) = {{ed^ + Q + dl'h^j) {{d^ + Q + dl)i,) 
=e(a,V)' + (QV')' + {dli,f 
+ 2{Qi,){dli,) 

+ (1 + e)(a2^)(QV;) + (1 + e)(52v>)(a^V)- 

In the last line, the contribution from d^ij) in Qij) gives a strictly positive term, so 
it can be dropped, and integration by parts can be applied to the remainder 

Ma2v)((Q + 5|)V) >^i{^',^p){lH') 

>^l\y^M^ + dt{ii{dt^){m - y • (M(5*v)(y5tV)). 

Finally, it remains to show that 

e|5,VP + m? + W + 2{Q^){dl^) ^e^l' + \{Q + dl)i;\' 

>e\dM' + mf- 

This follows from expanding the left-hand side using Q + = ^+ a'siv? Od"^, 
estimating the mixed term 2a^ sin^ ^(9j V)(^) by the Cauchy-Schwarz inequality, 
and using the fact that |ci| < e- □ 

It is sometimes useful to have an estimate which gives equality, except for some 
small error terms, instead of an inequality. In such cases, we relate Q + 3| to to 
observe that 

(£,v)(-cv) = {{ed't + m) ((^t + 

+ sin^ e {{d^i^){m + {{ed^ + ^)V)(a*V)) • 
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Thus, 



\{C,^P){C^) - {{edl + ^)^) ((a? + ^)^)| <a'\^M 



Expanding {{ed'^ + {{df + in the same way as in lemma 



2.4 



we find 



\{£^'ip){£^) — |V'l2,e + (time and angular derivatives)] <a^|7/;|2,;^, (2.4) 

where the time-derivatives terms are time derivatives of terms dominated by \dtip\ 1^24^1- 
Similarly, 

\{C^dt^f + {C^^ipf - IV'ls,^ + (time and angular derivatives)] <o^|V'l3.i' (2-5) 

where the time derivative terms are derivatives of quantities bounded by l^tTi?/)! |§2'iriV'| 
and where the angular derivatives are derivatives of smooth terms. 

2.4. Further notation. It is convenient to write the second-order symmetry op- 
erators with respect to coordinate partial derivatives 

Sa=-do.fiSfdi3. 

All other operator built from these, such as C, C^, and TZ, can be similarly expanded. 
For example. 



We use the notation 



/ = OirP) 



to mean that there is a constant, uniformly in a in some small interval of a values 
containing 0, such that Vr > r+, /(r) < Cr^. We introduce also the notation 



A 



g 



r-P 



to mean that there is a constant, uniformly in a in some small interval of a values 
containing 0, such that Vr > r+. 

This measures the decay rate at r+ and oo. If / is continuous, this is all the 
information that is required to bound the function. 

We use \x to denote the indicator function, which is identically one on X and 
zero elsewhere. 

We define a function to be smooth on a closed interval if it is smooth on the 
interior and if all the derivatives are continuous up to the boundary. 



3. The bounded-energy argument 

In this section, we construct a bounded energy by first constructing an almost 
conserved energy and then proving a Morawetz estimate to control the growth of 
the energy. 
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3.1. The blended energy. Recall from (1.8 1 that for \a\ sufficiently small, the 
vector field 

T^ = dt + x^nd^ 

is timelike in the exterior and Killing outside the region [r^ , + M] , since x is 
constant outside this region and decreases from one to zero inside this region. If we 
choose sufficiently large so that it corresponds to a larger value of r than any 
photon orbit for our initial choice of small |a|, this property will remain true for any 
subsequent decrease in the range of \a\ we allow. For specificity, we take — lOAf , 
which is beyond the range of photon orbits for any Kerr black hole. 

The vector field T^ becomes null on the horizon, so the associated energy degen- 
erates there. We compare this with the energy associated with T_l = (AI]/n)^/^nsj 
to make clear that the rate of degeneration with respect to the normal is roughtly 
(A/(r^ + a^))^/^. We also provide a coordinate expression which is useful for making 
estimates. The apparently singular contribution to the energy from (T ±ip)'^ is 
in fact vanishing, since the vector-field Tj^ vanishes on the bifurcation sphere at 
such a rate to exactly compensate for the factor of A~^, and then the form dr is 
degenerating at a rate of (A/(r^ -t- a^))^/^ near the bifurcation sphere. 

Lemma 3.1. There is a positive a such that for \a\ < a and any smooth function 
Tp, T^ is timelike and 



A 

(r2 + a2)2 



(T^V) + Hdri^) + Hdti^Y + 1®^^!'' (3.1) 



^ / T'^^rfV- (3.2) 

Furthermore, if is a solution of the wave equation Dt/i = 0, then 
1 



--AuH\drX\\dMdri^\- (3-3) 



Proof. We first expand 

A ' A 



( 'p'^ _L n^\^ An A/fr 



We now substitute dt = T^ — X'^ndij, and estimate the terms arising from the 
difference. We find 

Q''P{d^i^){dp^) =A{dAf - + °')' (T^V)' + g"''(5aV)(a,3V) + {d^^f 



- 2xc^H I {dttp){d^ip) 

+ ^ - (r' + a^u^h + AaMrxLOn) (9^) 



A 

^2 I „2\2 



--Aidri^r - ^° ^ (T^V)^ + g"''(5.V)(a/3V) + {d^^r 
+ aO{l, r-^){dti^){d^i;) + a^O{l, r-^){d^ijf. 
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Similarly 

+ a^O{l,r-^)dtij + aO{l, r-^)d^^. 
The t component of the momentum associated with T^^, is 



Since the asymptotics of the coefficients on the last line grow no faster than the 
coefficients of the terms in the first line, we can take a sufficiently small that the 
first line easily dominates the last, and we can conclude 

Since the term in Q has a bounded factor times a^, 



A 



i 

Since 



T_L =Tx + aO(A,r 



~3\ 



the difference between the two associated energies is 

which is easily dominated by V^^, and the two momenta are equivalent. Hence, 
their integrals are equivalent. 



We compute the divergence of the momentum using equation (2.11 and find it 
be given by 

which we estimate in absolute value. □ 
Recall that we defined higher-order energies by 

n 

£^T,,„+lM=Ei?Tj§.^], 
i=0 

where Si is the set of order-z symmetries from ( |1.7[ ). 
Corollary 3.2. If ijj is a solution of the wave equation Dip — 0, 

"TT-E'X n+l \drij\n\d4,ilj\nr'^drsmed0d(t), (3.4) 

at J^^ Js2 

where the norms on the right are defined in suhsection \2.!^ 

Proof. This follows from the previous lemma applied to the functions obtained 
from applying the symmetry operators to the solution ij: and then summing over 
the operators. □ 
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3.2. Set-up for radial vector fields and their fifth-order analogues. If z and 

w are smooth functions of r and the parameters M and a, then we can define the 
following single- and double-indexed quantities 

H- =—TZ- 
A ' 

n'- =dr (^n^ 



A 

AV2 ' 

These can be used to define a double-indexed family of vectors and scalars which 
we will use to prove a Morawetz estimate. 

Definition 3.3. The Morawetz vector fields and scalar functions are defined to be 

9a = 2 (^-^^ ) ~ 9a = 2 (^^^ ) - 9a" 

J'^ =zwn'^, T'^ =zwn'^^d^ , 

1% =\{drz)wn"^, qf, =]^{drz)wlZ'^'^d^. 

For simplicity, we introduce the following notation for the pair consisting of the 
previous vector field and function, 

A={A^,qf). 

Lemma 3.4. If tjj is a solution to the wave equation □■0 = 0, then the divergence 
of the momentum associated with these quantities is given by 



^V'^{Sai^){Sb^), (3.5) 



whe 



yob _y(a£b)^ 
V^=^{drAdrZ{drW'R.'^)), 



Proof. In the formula for the divergence of the momentum, (2.2 1, the terms in- 
volving Q"^{daSau){dpSbu) are refered to as the Lagrangian contributions. The 
Lagrangian contribution from (l/2)(9-yA''''''^'' 
contribution from the vector field. 



) in q^ exactly cancels the Lagrangian 
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The divergence of the momentum is given by 



^T^ (dr ) ) {d^Sa^){dpSb4^) 



1 

In the coefficient of the radial derivative terms, the part coming from the vector 
field can be rewritten as 



2 '" 7 V VAi/2_ 
Expanding using the definitions of z, w, and TZ' , we have 

A* 



+ ^{dpg"fd^{z{drWTZ'^^))C'-\Sa^j){SbjP). 

The expression TZ' was chosen so that it is exactly the derivative in the second term. 

Similarly, the quantity TZ" was chosen so that it is the derivative in the first term. 
Thus, the total bulk term is 

+ ^{drAdrZ{drWn'^'^))C'^\Sa1p){Sbtp). 

Since TZ"-C- is contracted against a quantity which is symmetric in ab, it is not 
necessary to distinguish between 'JZ"-C- and 'JZ"''-C-\ Substituting the definitions 
of A^, U^'^P, and gives the desired result. □ 

3.3. Rearrangements. We rearrange the terms related to U to get a strictly pos- 
itive contribution to the divergence. 

Lemma 3.5. If is a solution to the wave equation Dip — 0, then 

-da + =A'^{drSa^){drSbtP) 

where A, lA, and V are defined in lemma \374\ and 
We will refer to Ba-i as the first boundary term. 
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Proof. Starting from (3.5 1, it is only the second term on the right side that needs 
to be manipulated. First, we rearrange the derivative term to get 

The first term on the right can be rewritten in terms of Sc, which can be commuted 
with Sb, which in turn can be expanded in partial derivatives: 

-u'^c'-{Sai^){^^^iSf^pSbJJ) = - ^iU'^d^{Sa4^){s,Sb^) 

We can substitute this into the previous calculation, rearrange a derivative in the 
new expression, reindex, and use the symmetry of U— to conclude that 

- d^UjLl^C'-Sf{Sai^){dpS,i;)) 
+ dc.{fjLl^C'-Sf{Sai^){dpSbJ)) 

=^JU'^£"|'{^^Sai^){^0Sb1P) 

Applying the definition Ba-i gives the desired result. □ 

3.4. Choosing the weights. In this section, we choose the weights z and w to 
ensure the positivity of the highest order terms in the right-hand side of the estimate 
in the previous lemma, lemma [3751 



Definition 3.6. Given a positive value for the parameter eg2, we use the following 
weights to define the Morawetz vector field, 

A _(r2 + a2)4 

(r^ + a^)^' 3r2 — ' 

1 



This choice of weights generate a momentum which has a positive divergence, 
and for which this positive divergence dominates the square of third derivatives of 
tjj. The statement and proof of the following lemma make use of the norms given 
in subsection 12.31 



Remark 3.7. A simpler version of this argument can be run taking z = z\ and 
w = w\. These weights are chosen so that the coefficient of in TV and of d^dt 

in TZ" vanish respectively. Eliminating the &f term in TZ' is a natural first step, 
since, in essence, this is what has been done in all previous vector-field arguments 
in the Schwarzschild spacetime. Eliminating the d^dt term in TZ" is also natural, 
since this leaves only Q and d"^ terms, which have a natural ellipticity. Taking 
these choices gives a collection of vector fields for which the divergence is positive. 
However, these simpler choices generate a divergence which fails to dominate third 
derivatives which involve two t derivatives, and they generate an energy which is not 
dominated by the energy. We have further introduced the weights Z2 to obtain 
better control over the derivatives and W2 to "temper" the vector fields, so that 
the associated energy can be controlled by the energies associated with T^- 
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Lemma 3.8. There are positive constants a, eg2 , and C such that if \a\ < a and 
eg2 < eg2 and ip is a solution to the wave equation Dijj = then 

A2 ,^ 1 9Afr2 - 46M2r + 54M3 

>M 2J ^'-V'2,e,2 + R 4 

- c' ^2(^f+a2) ((« + 4?)|5.V'l2,i + esH^w^lLO 

-C^(a|V'||i + eo^HV'l2.aO- (3.6) 



+ aO{r'^)d^dt + a^Oir-^)Q + a^Oir-'')dl 
+ eo2a'0{r-')df + eo20{r-')Q + eo20{r-')dl 



and where the B'^.jj satisfy 

a a 

The angular components ofBA-ii are smooth. 

Remark 3.9. Here and through out he rest of the paper, we say that a vector field 
has smooth angular components if, for fixed r and t, the contraction of the vector 
field with any smooth 1-form on the sphere produces a smooth function. Because 
of the coordinate singularity in the {9, (j)) coordinates, this does not assert that the 
9 and 4> components of the vector field are smooth. The angular divergence of a 
smooth function has smooth angular components. In the applications of this lemma, 
it will always be the case that the regions of integration have boundary which are 
level sets of t or r, so that the boundary integral of Ba-ii does not contribute to he 
boundary integral. 

Proof. From lemma [STS] there are three terms to control, the A, and V terms. 

Step 1: The U term. The U term can be expanded using the definition in 
lemma 13.41 as 

so it is sufficient to calculate TV . With our choice of z and w, this is 
H' =- ea2(2(r - UI)r-^ + a^O{r-'^))dl 
+ aMO{r-*)d^dt 

- (2(r - 3M)r-^ + a^Oir-"^) + eQ20{r-'^))Q 

- (2(r - 3M)r-4 + a^Oir-") + e920{r-'))dl. 
Step 2: The A term. The A term is 

— I — 
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With our choices of z and w, we find: 



- aMO{r-^)d^dt 

+ M{r-^ + a^O{r'^) + e920{r-^))Q 
+ M{r'^ + a^O{r-^) + eQ20{r-^))dl. 



We are interested in this because the operator TZ" is very close to C^^^ in the 



sense that 



2 '~'<^o2 



--aO{r-^)\Tldr^l^\ 

+ aO{r-^)\Tldr^\ 
+ aO{r-^)\Tldrij\ 

+ ea20(r-3)|0?9^VI + eg2aO{r-^)\sm'^ 9dfdr^l^\ 
--aO{r-^)\drij\2,i + el20{r-^)\drij\2,i + eQ20{r-^)\drij\2.a^. 



Since C and >Ce ^ commute with functions of r, we can apply lemma 2.4 to dri^, 

I I 

to get 



M 



j.2(^j.2 _|_ g2^ 



{C,.dri^){Cdri^) >M 



A2 

^2(^2 + a2)l^-V'lL,. 

time and angular derivatives. 



The time and angular derivatives are exactly those coming from lemma |2.4[ The 
terms from the angular derivatives are smooth, and the terms from the time deriv- 
ative are of the form 



M- 



dtmdr^Ki4^+ a sin^ edl)dr^k) 



Thus, we only need to control contributions from these terms when they appear as 
boundary terms on hypersurfaces of constant t. They are controlled by 



M- 



■^(r^ + a?') 



\dtdr^k\\{/f^+ sin' edl)dri})\ <M 



^2(j,2 _|_ q2^ 



a 

(3.7) 



Thus, 

A^{Sadri'){Sbdr^) >M 



j,2^^2 _|_ q2^ 
A2 



C ,A a\drMil+^0^\drMia^+eUdr^\2,l) 

+ * ' * 



r2(r^ _l_ Q,2^ 
time and angular derivatives 



with the time and angular derivatives satisfying the bound given in the statement 
of this lemma. 
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Step 3: The V term. By direct computation, the V term is given by 

V^Sa^^drAdrqlSa 



+ aO{r-^)d^dt 



1 

+ {a + es.)0{r-^)§2, 



where we have used 0{r '*)§2 to denote terms of the form 0{r ^)Sa with Sa G §2- 
Applying the estimate in (2.4), we find 

V^{Sai^){Sbi^) =^(9Mr-2 - 46AfV-3 + 54AfV-4)(/:^^)(/:^) 
+ 0(r-2)(a|^|2 , + eg.|^|2 J 

>^(9Mr"2 - 46AfV"3 + 54AfV"^)|V'|2,eg2 

+ 0(r-2)(a|V'|^,i + ea2|^l2.a) 
+ time and angular derivatives. 



Again, the time and angular derivatives come from the application of lemma 2.4 



so that the terms from the angular derivatives are smooth, and the terms from the 
time derivative give a contribution of the form 

^\dtM4^+a'sm^9dl)^\ <^|9t^|^|5aV'l- (3.8) 

a 

The time and angular derivative terms arising in this step and the previous one 
are combined into Ba-u and are controlled by and (3.7l-(3.8l. □ 

Lemma 3.10 (Controlling the boundary terms). If ip is sufficiently smooth, sat- 
isfies Oip = 0, and has initial data which decays sufficiently rapidly at infinity, 
then 

\B%j\ <C(7'tJS2V]* + |§2^n 

\bXjj\ <c{VTjdt^Y + VTjS2i^Y + \dttp\^ + l§2^n, 

lim Va =0 

r — »-r_|_ 

lim Va =0 

r — >oo 



Here, by "sufficiently rapidly", we mean that the T^ energy is convergent. In par- 
ticular, lim^oo "0 — 0; hmr_>oo rdrtp = 0, lim^— >oo rdtip = 0, and the same estimates 
hold for Slip and S2'4'- 
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Proof. By direct computation, 

= (- ((r^ + a2)2 _ aZ^sin^e) (8^2^) - 2aMr(9^§2^)) (0(l)(9r§2^) + 0(r-2)(§2^)) 
=0(r-i)((§2^)(Ti§2V)) 
n =0((A/r2)2,r2)(a,S2^)(9.§2V') 

+ 0(A/r2, r)(a,§2V')(S2^) + 0{A/r\ 1)(§2V')(S2^) 

Here we have used S21P to denote a term which is bounded in absolute value by 
1 82'!/' I- First, the t component of the momentum can be dominated by 

+ 0((A/r2), r2)(a,§2^)' + 0(A/r2, 1)(§2^)^ 



From the expression for the energy, (3.11, the second term is dominated by [§2'/'] • 
The third is clearly dominated by ip'^. This leaves the first. Since the coefficients 
in the first term are rational, and the ratio of the coefficients of the dt^ and d^ip 
term in the first expression is 

{r^ + aY -a^^sin^e {rl + a^^ 



2aMr 2aMr+ 



+ 0(A/r^r^) 



a 

=c.^i + 0(A/r2,r4), 

the first term can be estimated by 

c|-(atS2V' + i^Hd^S2^f + 0{r''){dS2^f + 0{l){d^S2^f 



Thus, the expression for the energy in estimate (3.1 1 shows that this is bounded by 
Vi^ [§2^]- This controls the t component of the momentum. 
The t-componcnt of the first boundary term is controlled by 

B%i =0(r-i)((S2V^)(Ti§2V')) 

<0{r^){dt^2i^f + 0(l)|Oi§2V'|' + 0(1)(S2V)' 

<7'^J§2^] + |§2V'r 

The t component of the second boundary term was partially estimated in the 
statement of lemma 13.81 so that 



+ 0{r-^)\dtn^2'>P\ 
<0((A/r2)2, l)\drdM^ + 0((A/r2)2, l)|a,§2^P 

+ |5t^P + |§2^P 

<C(7'^J5tV] + ^,[§2^] + \dM^ + |§2V^n. 

The limits at r+ and 00 are easily evaluated. The radial component of the 
momentum consists of bounded functions times a power of A, so they vanish at 
r+. From the asymptotic behaviour of the coefficients T— and g^, which are 
dictated by the asymptotics of z and w, it is sufficient that 0{r'^){dtStp){drSip) + 
0{r){drSilj){Sip) + 0(1)(S'V') converges to zero, at least on a sequence. This is 
a necessary condition for the energy to be bounded, so it imposes no additional 
restriction. □ 
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Note that it is not necessary to estimate the hmits of the radial components 
boundary terms, S^.j and iB^.u, since these components are identically zero. 



3.5. The Hardy estimate. In (3.6l, the coefficient of IV'li e 2 positive except 
in a compact range of r values. The purpose of this section is to prove a Hardy 
estimate which allows us to get a globally positive coefficient for l-iAli.e 2 using 

the positivity of the term involving ISr-V'li e 2 ■ '^^'^ proof is a bit technical, and the 
reader can omit it on a first reading, since the proof is independent of the rest of 
the Morawetz estimate. 

Lemma 3.11. There exist positive a and euardy such that if \a\ < d, then for any 
smooth function (p on [r-|_,cx3) x S"^ which is bounded on [r+,00), 

A2 ^I9r^ - 46Mr + bUP 

r+ yr2(r2 + a^) 6 J 

roc ^2 ^ 

>euaray ^^^^^^{dA? + ^^'dr. (3.9) 

Proof. The proof consists of several parts. The early parts of this proof follow the 
method of |5]. First, we will demonstrate that it is sufficient to find a positive 
solution to an associated ODE (ordinary differential equation). Second, we rewrite 
the estimate and ODE in terms of a new function, ip. Third, we will construct an 
explicit solution for the new ODE when a = and euardy = 0. Fourth, we will argue 
that the construction of the explicit solution can be perturbed to cover nonzero a 
and CHardy, which will give a perturbed estimate for ip. Fifth, we will show that this 
gives the estimate for the original function (jj. Finally, we will show that boundary 
conditions for the ODE do not place restrictions on the function (f>. 

Step 1: Find a positive solution to the associated ODE. We wish to 
show that if the ODE 

-drAdrU + Vu^O, 

has a smooth, positive solution u on [ro,cxD], then for any smooth function </> on 
[rg, 00], there is the estimate 



00 



A{dr(l3f + V(j)^dr>0, (3.10) 



'•0 



as long as 



cjj^A^ (3.11) 
u 

vanishes at tq and 00. Recall that we define a function to be smooth on a closed 
interval if it is smooth on the interior and all derivatives have a limit at the boundary. 

Since u is positive, for any smooth 0, we can define / = (p/u. From integration 
by parts. 



00 /-oo 

2 , T^j,2j„ r /!„, „, -fMoo _ / „,f2/ 

ro 

ro •'ro 



A{dr(l)f + V(p^dr - [Aufidruf)]^^ = / uf^{-drAdrU + Vu)dr 



00 



u^A{drffdr 



~ [u'Af{drf)]r,. 

Since u satisfies the ODE —drAdrU+Vu — 0, the first term on the right is positive. 
Cancelling the boundary terms on the right from those on the left leaves the estimate 

/>oo 

A{dr(t>f + V(t)'^dr = / u^A{drffdr + [f^ Au{dru)\^^. 

ro Jro 
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The boundary term vanishes under condition (3.11 1, and the integrand on the right 
is non-negative, since = fu. Therefore, 



A{dr 



Vcffdx >0. 



Step 2: Simplify the estimate to eliminate one of the coefficients. We 

will first introduce the notation 

We will consider the function 

Since A^/"^ is smooth on [r+,cx)) and vanishes linearly at r^, the new function {p is 
also smooth and vanishes at least linearly at r+. Its derivative satisfies 



1 



1 drA 

Therefore, the left-hand side of ( 3.10[ ) is given by 



drA 
~A~ 



(p{drip) 



I jdrAf V 

A A^ A 



Y_ 1 dlA 1 



Lp dr 



93^dr. 



1 drA 



If the original function </) is bounded, then the boundary term in this equality 
vanishes. If the integral satisfies 



V \dlA 1 (drAf 



^A 2 4 A2 

then by multiplying this estimate by 1 — enardy.a 
'V 18^ A l{drA)^ 



tp'^dr >eHardy,i 



1 



(fi dr 



(dr^y 



A 2 4 A2 



(f dr 



> 



V 1 A 1 {drAf 

eHardy,3 I ^ + 2 ^ " 4 



(1 - eHaidy,3)eHardy,2^;^¥' ) dr. 



By taking enardy.a sufficiently small, we can conclude inequality (3.9) holds. 

Step 3: Construction of the explicit solution for a = and enardy = 0. 
Following the arguments in the first section, we could prove the desired estimate 
(for a — and enardy = 0) by finding a positive solution to 

-drAdrU + Vu = (3.12) 

with 

(7-2 - 2Afr)2 



A 



V = 



1 9r2 - 46Mr + 54M^ 



6 r4 

on the interval [2M,oo). However, by using the argument in the previous section, 
it is easier to use the transformed function 



X =r - 2M, 



(3.13) 
(3.14) 
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and to solve the ODE ( |3.12 l 

-dlv + Wv^Q, (3.15) 

V ldlA_\{d^ 
A 2 A 4 ^2 
_9x^ -UMx-2M'^ 

6a;2(a; + 2Af)2 ^ ' ' 

on the interval x G [0, oo). 

We first note the following properties of hypergeometric functions [TJ [12]. The 
hypergeometric function is typically written with parameters F{a, b; c; z). It should 
be clear in all cases whether a refers to the first parameter of the hypergeomet- 
ric function or to the angular momentum parameter of the Kerr spacetime. The 
hypergeometric function F(a, b; c; z) has the following integral representation for 
a < < b < c and z ^ [1, oo) 

r(c) /"^.b-lj-i j-\c-b-l/ 



F{a,b;c;z)^ ^ / t'-\l - 1)^-'-^! ~ tzy^dt. (3.17) 

r(a)r(6) Jo 

It is not obvious from this representation, but it is true, that F is symmetric in its 
first two arguments, F{a, b; c; z) = F(b, a; c; z). There are a vast number of further 
relations. The hypergeometric differential equation is 

d'^w dvj 
z{l - z)—Y + [c- {a + b+l)z] abw^ 0. (3.18) 

A pair of solutions to this equation is 

F{a,b; c; z), 

z^-^Fia - c + 1,6 - c + 1;2 - c;z). 

Returning to the ODE arising from the Hardy estimate, we introduce the further 
substitution 

v=x°'ix + dfv. (3.19) 

The ODE now becomes 

v" = (a(a - l)x'^-\x + 2f + 2a[3x°'~\x + df-^ + /3(/3 - l)x°'{x + df-"^) v 
+ 2 {ax^'-^ix + df + Px'^ix + df-^) v' 
+ x°'{x + dfi", 
= -v" + Wv 

{x + df-^P, (3.20) 



Q-2 



=X 

P =x^{x + dYi" 

- 2x{x + d){{a + [3)x + ad)v 

+ ^ - {a{a - l)(x + df + 2al3x{x + d) + P{(3 - l)^^) 

9X^ -UMX -2M^ 2/ n2 \~ 



' 6.2(. + 2M)2 + (3-21) 

The prefactor of x°'-'^{x + df-"^ in the ODE can be ignored. If we choose 

d ^2M, (3.22) 

then the rational function in the last term on the right reduces to a polynomial. 

The coefficient of v" is x'^{x + d)'^, of v' is x{x + d) times a linear function, and 
of w a quadratic. If we choose the parameters a and /? so that the coefficient of 
■D is a constant multiple of x(x + d), then an over all factor of x{x + d) can be 
dropped, leaving the coffecients of v", v' , and v as x{x + d), a linear function, and 
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a constant respectively. The substitution z ~ —x/d, then transforms the equation 
to the hypergeometric differential equation. Our goal is to show such choices of a 
and (3 can be made. 

It is now merely a matter of checking by direct calculation that this can be done. 
The coefficient of v is 

- a{a - + 2xd + d^) - 2af3{x'^ + dx) - (3{f3 - l)x^ 

+ (3/2)a;2 - (17/3)Mx - (3.23) 

In this coefficient, we set the constant order term to zero 

-a{a - 1)^2 - m2/3 =0, 

1 ^ V6 
"=2^^- 

Fortunately, the term aP{x'^ + dx) is already a multiple of x"^ + dx, so we may ignore 
it when trying to get the coefficient of u to be a multiple of x'^ + dx. We set the 
ratios of the remaining coefficients of x'^ and of x to be d: 

d ((3/2) - a(a - 1) -/?(/?- 1)) - 2da{a - 1) - (17/3)M. 

We can substitute —a{a — 1) = 1/12 to get 

2 ((3/2) + (1/12) - /3(/3 - 1)) -(1/3) - (17/3), 

^2 2 

The four choices of sign provide four choices of simplified equations to study. For 
simplicity, we will consider only the equation arising from taking the + sign in a 
and the — sign in /3j^ 

We are left with the differential equation for v 

x{x + 2)v" - 2((1 + \/6/6 + 3V2/2)x + 1 + V6/3)v' 

+(19/6 - 3V2/2 + V6/6 - VS)v = 0, 

Making the substitutions z — —x/d and ip{z) = v{x) gives 

z(l - z)'iP" + ((1 + \/6/3) - (2 - 3^/2 + %/6/3)z) ^P' 

+ (-19/6 + 3\/2/2 + \/3-\/6/6) 7/^ = 0. (3.24) 

Thus we have a hypergeometric differential equation, with solution v = F{a, b; c, —x/d). 
We can immediately read off some quantites in terms of the hypergeometric param- 
eters 

c=l + \/6/3, (3.25) 
-a-6-l=-2 + 3V2 - 76/3, 

-ah = - 19/6 + 3\/2/2 + \/3 - \/6/6. 
We can now solve for the remaining two parameters 

{aM=\-l^2+^±\V7. (3.26) 
We will make the choice a < 6 so that 

a < -2.5 < < .1 < 5 < .2 < 1.8 < c. 

In particular 

a < < b < c. 



^This choice simplifies some expressions in the rest of this argument. 
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Thus, the integral representation (|3.17p holds. Multiplying by (-l)"'T(c)/(r(a)r(&)), 



we find that '4'{z) is positive when z < 0. This means that v is positive when a; > 0, 
V is also positive when x > 0, and u is positive when r > 2M. 

Step 4: The perturbed estimate for v. In this step, we will prove that 
there are < aHardy,4 and < eHardy,4 such that for \a\ < aHardy,4 and all suitable 



for 



+ Wip'^dx > 0, 



^-^ _ 34^^ _ (M + x)2 

W — :r^77 eHardy,4^ 



6x2(x + rf)2 

d =r I — r_ 



x2(x + dy 



This potential is of the form 

.T, ^ Cix2 + C2X + C3 
C74x2(x + d)2 ' 

with the coefficients Ci, . . . , C4, and d perturbed from their original values in equa- 



(3.27) 



tion (3.161 



From the argument in step 1, it is sufficient to find a positive solution to the 
associated ODE ( |3.15| , 

with the perturbed potential W. The analysis in step 3 found an explicit, positive 
solution for x € [0, 00) for the parameter values dictated by the potential in equation 
(3.161. This step shows that the previous analysis also applies when the coefficients 



are perturbed. 

The previous analysis began by making the definition of v in equation (3.19), in 
terms of the parameters a and p. The analysis then proceded by choosing values 
for a and (3 by solving quadratic equations coming from the coefficient in formula 
(|3.23|, which lead to the new ODE (|3.24|). This ODE could be solved explicitely 



in terms of a hypergeometric function by solving linear and quadratic equations for 
the non-zero quantities a, b, and c. Since the coefficients in formula (3.23) depend 
continuously on the parameters Ci, C2, C3, C4, and d in the potential; since the 
coefficients in the ODE (3.24) depend continuously on a, /3, and the coefficients 



in the potential; since all the quadratic equations involved had distinct, real roots; 
and since solutions to linear and quadratic equations depend continuously on the 



coefficients; it follows that positive solutions to the ODEs (3.15) and (3.24) can be 
found explicitely in terms of hypergeometric functions with parameters a, b, and 
c depending continuously on the parameters in W, at least when those parameter 
values are sufficiently close to the values given in equation (3.16). Similarly, when 



the perturbation of the parameter values in the potential W is sufficiently small, 
then the hypergeometric parameters maintain their order a < < b < c. This gives 
the existence of positive aHardy,4 and eHardy.4 which give the desired estimate for 
this step. 

Step 5: The perturbed estimate for the original function 0. In the 

previous step, a particular type of perturbation of the potential was considered. 
In this step, we show that such perturbations are sufficient to control the type of 
perturbation appearing in our probelm. 
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From the argument in step 2, we wish to prove that there exist < a and 
< eHaidy.2 such that for < |a| < a and suitable ip 



with 



2 4 



dr >eHardy,i 



1 1 

j,2 



ip'^dr, (3.28) 



V 



1 9r2 - 46Mr 
'6 ^ 



54M2 



With the benefit of foresight, we introduce a new rotation parameter 

a =Af - - a?. 

When a is treated as a function of |a| with M fixed, this is a continuous, increasing 
function on the interval [0, A/], which maps the interval [0, M] to [0, A/]. In addition, 
since the quantities which appear in our estimate (such as A and r'^ + a?) only have a 
quadratic dependence on a, and since a? can be solved for as a quadratic expression 
in a, it follows that the quantities A and V are rational functions in (r, Af, a). 

The new radial coordinate, analagous to the one defined in (3.141, is now defined 
to be 

X =r - r+ = r - {2M - 2d). 

Since r can be solved for linearly in terms of {x, A/, a), the quantities A and V are 
rational functions in (r, Af, a) . 
The quantity 

V Id^A _ 1 (drA)'^ 
2~A ~ 



A' 2 A 4 yl2 

is rational in (x, Af, a); has degree, with respect to a;, two lower in the numerator 
than in the denominator; has singularities in x e [— d, oo) only at x € {0, — d} for 
fixed M and a; these are of order at most two; and, for sufficiently small a, has no 
singularities in a for fixed a; > and M . Thus, we may expand it as 

A2 Qo + aQ>' 

where the functions Pq and Qo polynomials in (x, Af), the functions P> and Q> 
are polynomials in (x, Af, a), and Qo and Qo + aQ> have no roots in x G [— d, oo). 
Since Pq/Qo is determined explicitly by equation ( 3.16[ ), it follows that 

^ ~a PyQo-PoQ> 
A^Qo A^Qo{Qo + aQ>) 



W 



must decay like r as r ^ oo for fixed d and Af and has no singularities in [—d, oo). 
Since this is a rational function, there is a constant C such that 



W 



<dC 



(Af + a;)^ 



A2 



J^Pa^ 
A2 Qo 

Thus, there are sufficiently small d and eHardy,2 such that for < \a\ < d 

1 



W- 



eHardy,2^ > W, 



with as in equation (3.27) The smallness of d and eHardy,2 is determined by the 
smallness of aHardy,4 and eHardy,4- These then give d and enardy for which the desired 
estimate holds. 
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Step 6: Controlling the boundary terms. Since the argument from step 1 
was applied to the function Lp, the boundary condition which must be imposed for 
this argument to hold is that 

V — 

V 

vanishes at r_|_ and at co. Since the positive solution to the ODE is given by 
v{r) =a;°(x + dfi = x°'{x + dfF{a, b; c; -z/d) 
=(r - r+nr - r^fP [a, b; c; -^^^^ . 

and the hypergeometric function is analytic (in its fourth argument) near zero, the 
ratio drv/v will diverge at most inverse linearly at r Thus, it is sufficient 

that if vanish linearly at r = r+. Since (/? = A/(r(r^ + a?Y/'^)(j), it is sufficient that 
(j) be smooth near r+. 

To show the decay near oo, we first note that from the form of the potential W in 
the ODE, the solution v{r) will behave like a polynomial as r ^ oo, so that drv/v 
will decay like a constant times l/r. Thus, it is sufficient that Lp remains bounded 
at infinity. Thus, to get decay at both r+ and oo, it is sufficient that <j) be smooth 
and bounded on [r^-jOo). □ 

3.6. Integrating the Morawetz estimate. 

Lemma 3.12. There are positive constants a, f, Ci, and C2 such that, for all 
\a\ < a and all smooth ip solving the wave equation Dip = 0, the estimate 

Ci{Et, [§2^1(^2) + Et, [SiV'](T2) + St, [§2^'](Ti) + Et^ [SiV'KTi)) 

fT2 /-oo /. /a2 



> [ ^ r [ (^)\drMh + r''\^\ll + ^r^3M-mlld'fM (3.29) 

JTi Jr+ JS^ \ ^ J 1^ 

/'T2 roo f- 1 

-^C2 / / / lr^3M-|^|2rf'Ai, (3.30) 



Ti Jr+ 



holds, where lr^3Af is identically one, except in an open neighbourhood of width f 
about the values of r for which there are orbiting geodesies. In this neighbourhood, 
we take lr-^3M to be indentically zero. 



Proof. We integrate the result of lemma 3.8 over the coordinate slab {t, r, 6, <j)) S 
[ri,T2] X (r+,oo) X S'^, from which we get the integral of the right-hand side of 



estimate (3.6). From the Hardy estimate (3.9 1, the integral of the first two terms 



on the right-hand side of (3.6 1 dominates an absolute constant times 



By taking \a\ < eg2 < 1, these terms will also dominate the fourth and fifth terms, 
with a constant factor left over. Since eg2 can be chosen independently of a, the 
norms \1p\2 e 2 can be replaced by \ip\2,i at the price of a fixed constant. The same 
is true for the norms of dril^. 

The only term which we still need to estimate is the third, 

T2 foa n (- 2 I 2\4 

, i Is. 4!:(3^3i^^"'(^"^>)(^^^>)dV. 
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The integrand can be estimated by 

>Clr^3Mr-^\riC,.,i;\^ 

+ l,.^3MO(r-5)eg2|Ti/:,2V'P 

+ lr^3MO(r-')(|a| + e|2)|TiS2V|'- 



The first term in this expression can be bounded from below by inequahty (2.5 1, so 
that 

|Ti£e.2'/'P Jil'/'la.e 2 + ^^^IV'lli + (time and angular derivatives). 

To control the remaining terms, we note that 

ea.|Ti/:,.VP + (|a| + 6|.)|Ti§2V'l' lii^a^a' + \a\ + e|2)|9?^P + {e^^y + \a\ + e^^Mm 

+ (eg. + \a\ + el,)\dt/H? + (^Of + l«l + 4?)l Wl 

If we impose the conditions that \a\ < ^ 1, then these terms are controlled by 
I^IIj ^. This norm can be estimated by (2.5 1 after commuting angular derivatives 
through Cf^, 2 . This calculation gives: 

IV'ls e 2 — I'Ce. 2'^i^P + time and angular derivatives 

<|£e 2'^iV'P + a^|V'l2 + time and angular derivatives, 



2 

2 



so that 

(r^ + a^)^ 
2r(3r2 — a^J 

+ lr^3A/0(r^"'^)(time and angular derivatives). 

Thus, if we fix a sufficiently small eg2 , and then require |a| < 6^2 , we have control 

of |w|3 with the weight ^r^3Mi^~^ ■ 

The time derivative generated in this part of the argument is 

at(1.^3AfO(r-i)(atTiV')(^Ti^)), 

where Ti is the set defined in section |2T4| to contain dt and the rotations around the 
coordinate axes. Thus, the contribution of this time derivative on the boundary of 
the region of integration is bounded by [SiV'] + Pt [§2^']- 

We must now control the integral of the momentum and the boundary terms 
over the boundary of the slab. All the angular derivate terms vanish, since has 
no boundary. Similarly, the boundary contributions along r — and r 00 are 
zero by lemma 3.10 (Geometrically, one would expect this, since r = r+ is actu- 
ally a two-dimensional surface, the bifurcation sphere, and not a three-dimensional 
hypersurface, so it should not contribute any boundary terms.) 

We are left to control the integral of the momentum and the boundary terms 
over the hypersurfaces t = Ti and t = T2. From lemma p?.10[ these are controlled, 
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at fixed t, by 



n + 6^,i + 6ViidV 



'St 

<^^Tj§2V'](t) + i^TjatV'](t). 

In this, we have used the 1-dimensional Hardy estimate: 



\^p\^dx< / dx (3.31) 

with X — r — r+ to control /^^ IV'pd^/z by ET^{t). □ 

We are unable to make use of the previous lemma since it controls only third 
derivatives, but the boundary terms involve both the second- and third-order ener- 
gies. (Certain second-derivative terms are controlled, but these are not the impor- 
tant ones.) In the following lemma, we control the lower-order derivatives. This is 
the only place where we use any separability or spectral information. 

Lemma 3.13. There are positive constants a, Cgp f, and C such that for all \a\ < a 
and all smooth ip solving the wave equation □'0 = 0, the estimate 

C{Et^.M{T2) + ET^^MiTi)) (3.32) 

- IJ r IsX (^) ^ ^ ^'""'''^r ^'^'"^'^ + '^^'^) ) '^M- 

holds, where lr^3M is identically one, except in an open neighbourhood of with f 
about the values of r for which there are orbiting geodesies. 



Proof. To control the higher energy, by corollary |3.2[ we must control the weighted 
integral of the first, second, and third derivatives. The Morawetz estimate, lemma 



3.12 controls the third derivatives, so we only need to control the lower-order 
derivatives from corollary |3.2| and from lemma |3.12| itself. 

Since commutes with the d'Alembertian, it is possible to separate out the zero 
eigenmode and write ip in the form 

■0 = ■0i.=O + ■0L,5^O- 



From the Morawetz estimate, (3.6 1, we control the integral of the third derivatives. 
If ip is real valued, then ipL^^o and 4'l^=^o are also real valued. 

Since the eigenvalue, L^, satisfies \Lz\ > 1, for tpL^^^o, the integral over the 
sphere of the third derivatives controls all lower derivatives, i.e. the homogeneous 
norm controls the inhomogeneous norm: for n = 2 or n = 3 (for which we have 
defined the homogeneous norms), 

/ IV'L.^ol^dV <C f I^L.^ol^.ldV- 



By taking \a\ small, this means that the second-order norms in the Morawetz esti- 
mate |3.12| can be dominated with a small loss by the third-order norms. 
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It remains to control the lower derivatives of ipL^=Q- To do this, we prove a 
Morawetz estimate using a classical, first-order vector-field. Consider the momen- 
tum associated with 

A =zwfdr, 



q'a {drz)wf, 



Here / plays the role of TZ' . Using the same sort of calculations as before, we can 
obtain the analogue of (3.51 

with 



Ai/2 



V ^^drAdrZ{drWf). 

Taking the same choices of z and w as before, we find 



A = 



1 A2 



2 r'^ + \r 



+ \a\Oir--') + +eo20{r-') 



A 

' (r2 -I- a^)^ 
1 9Mr2 - 46APr + bU'P 
6 



1 - 2e«2 



A 

(7-2 + a'^y 



+ {a + eof)Oir-^). 



Note that there are no dtdcf, or 9^ arising from TZ' when acting on 'ipL^=o- From the 
Hardy estimate, (3.9 1, it follows that 



Aidr^PL.^of + V\^JjL^=o\' 



>. 



\dri^L.=0\' + \\i^L.=0\'. 



Thus, 



> 



A2 



9*V ~r2(r2 -t-a2) 



In analogy with the previous results in lemma 3.10 there is a constant and an 
upper bound on a such that 



|i?(A,,.)[^L.=o]l<C^i?TjV'L.=o]. 
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Thus, we have the Morawetz estimate 

-L 



which controls the first derivatives. 

Second derivatives have either two, one, or zero time derivatives. 



(3.33) 



(3.34) 



Any term 



containing one or more time derivative can be controlled by applying (3.34) to 
dt4^L^=o- Terms which contain no time derivatives contain at least one angular 
derivative. Thus, by integrating by parts in the angular derivative and applying 
the Cauchy-Schwarz inequality, we can control such terms by the first and third 
derivatives, which are controlled by (3.34) and (3.6). This argument applies to 
lower derivatives arising from either corollary |3.2| and from lemma |3.12[ In the 
case of the lower-derivatives which must be estimated for the Morawetz estimate 
in lemma 3.12 we can use the smallness of \a\ to obtain the desired estimate for 
ipL^^a in this lemma. 

By combining the results for iJjl^=o ^^nd il'L,^ai we have the desired result. □ 

3.7. Closing the argument. We are now able to show that the energy associated 
with T;^ is uniformly bounded by its value on the initial hypersurface. When a = 0, 
the energy is conserved. When a 7^ 0, the energy is no longer conserved, but the 
factor by which it can change vanishes linearly in \a\. 

Theorem 3.14. There are positive constants a and C such that if \a\ < a and ip 
is a solution to the wave equation Dip = 0, then for all t2 > ti > 0: 

ET^^sMh) < (l + C|a|)^T„3M(ii). 

Proof. By lemma [3?2| 

<\a\C 



[ti,t2]x(r+,co)xS2 



Isuppx' (l'5rV'l2 



1^13) dV 



By the Morawetz estimate, lemma (3.13 1, for sufficiently small a, there is a constant 
C" such that the integral of the third derivatives is controlled by the energies. Thus, 

ET,,Mit2) - ET^^Mih) <\a\C' {ET^,3mt2) + ST,,3[V'](tl)) ■ 
Thus, for sufficiently small a, 

(1 - |a|C')i?T„3M(i2) <(1 + \a\C')ET,Ami), 



ET,,3mt2) < 



1 



\a\C' 



1 - \a\C' 



Since, for sufficiently small \a\, the rational function (l + |a|C")/(l — |a|C") is bounded 
above by 1 + C|a| for some C, the desired result holds. □ 



4. Decay estimates for the local energy 

In this section, we prove decay using the K vector field. To prove decay, we find it 
convenient to work in a different set of coordinates and to work with a transformed 
field, ■(/'• This allows us to write the wave equation as the Euler-Lagrange equation 
for a different action. A single equation is in general the Euler-Lagrange equation 
for more than one Lagrangian. The advantage of the action chosen here is that it 
allows us to control terms involving 1-011 with the energy. 
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4.1. Reformulation of the problem. Recall that the r» coordinate is defined by 

dr _ A 

and by the choice of initial condition that r» = at r = 3M. In the Schwarzschild 
case, this choice of initial condition is convenient since it fixes the origin of the 
coordinates at the photon orbits, and, in the Kerr case, no other choice seems more 
convenient. In this section, we will work with the (t, r*, 6, 0) coordinates and treat 
r as a function of . All indices will refer to these coordinates. 

Lemma 4.1. Let 

'ip —\/r'^ + a^ip ■ 
The wave equation Oip — is equivalent to 

Ai^'a, (m^"") 9/3V^ - V^V^ -0, (4.1) 

where 

V ^ , , ^ .,A '2Mr^ + r^a^ - AMra^ + a^). 
(r2 + a2)4^ ^ 



Here, it is understood that TZ is the operator defined in (1.6 I. 

Proof. Direct computation. □ 

The matrix G^f^ is obtained from multiplying the inverse metric g°'^ , evaluated in 
the {t, r^,,9, (f) coordinate system by the function SA(r^ + a^)~^. The transformed 
wave equation in (4.1 1 can be written as the Eulcr-Lagrange equation for 

S = CdL^ii^, 



£=1 (g"'3(a„^)(9/5^) + I/^2^ , fi^smt 



Recall that d^/i* was defined at the start of section [Sj Section 2.1 defines the 
canonical energy-momentum tensor, the momentum vector, and the energy for this 
Lagrangian. We denote these 

The related bilinear terms are defined similarly. 
In the r^ coordinates, we have 

ga0 ^ _ 7v-2t'^T^ + + h"'^, (4.2) 

with 

N-'^ ^l-a^sin^eVQ, (4.3) 

, T> f, 2 ■ 2fl'^^ + 2Mr + a2cos26'\ 1 ,^ ^, 

+ h - a' sm^ j — 2^^?^0- (4.4) 
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4.2. Energies in the reformulation. 

Lemma 4.2 (Dominant energy condition for T). There are positive constant a and 
C, such that, if \a\ < a, and if ^ and Y are future-directed, non-spacelike vectors 
(with respect to dt and g), then 

-f-^ 0g^^X.^Y^ > (4.5) 

In particular, 

CV^r^ > iT^i,f + (drjj)' + Vgiy^p + V\i^\\ (4.6) 
Proof. The energy-momentum tensor can be written as 

-fii^r^ =^"^(a^^)(<9^^) - ^-s'^pG'^idMds^') - Is'^pv^'- (4.7) 



AT 1 



Thus, the reformulated energy-momentum tensor can be written as the sum of two 
terms. The first satisfies the dominant energy condition because it is a positive 
multiple of the standard energy- momentum tensor evaluated on ip. The second 
term is a negative factor multiplied by the metric, so it also satisfies the dominant 
energy condition. Thus, Tlij}]"^ pg^^ is the sum of two terms which each satisfy the 
dominant energy condition, so it satisfies the dominant energy condition. 
We now turn to the second part of the lemma. From equation (4.7 1, 

t"0T^(di)„ = - (T^^)^*^(9^V^) -t- \g'<\d^^){ds^) H- \v^^ 

From the fact that T^(9„ ^ = T_l + aO{{I\/r'^),r-^)d^, and the 

expansion of Q in (4.2 1 and (4.4 1, 

r"/3T^(di)„ =7V-2(T^^)(Ti^) + aO{{^/r^),r-^){T:^ij){d^i,) 

+ i (-iV-2(T^^)2 + (9,^)2 + h'^f^{dj){dpi,)) + ^-Vi^\ 

Since we assume that a is small, and since 0{{A/r'^)'^ ,r^^) decays more rapidly 
than the product of 7V"2 ^ i and 0{{A/r'^),r~'^), which is the decay rate of the 
leading-order piece of h, the term aO{{A/r'^) ,r^'^){T ±'ilj){dci,ip) is controlled by the 
{l/2){T^'ifj)^-\-h"'^{dai^){dp^)) terms with only a small loss. This proves the desired 
result. □ 

Lemma 4.3 (Relation between energy for u and ip). There are positive constants 
a and C, such that, if \a\ < a, and if 4' = (''^ + a'^Y^'^tp^ then 

(1) (Global identity) 

ETj^P]{t)^ETM{t)- (4.8) 

(2) (Local estimate) Given ti < r2 and 7'*]^<— l<l<r*2 such that r» = r^^i 
corresponds to r — ri and = corresponds to r = r2. 
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Proof. This argument follows by substituting -0 = (r^ + a^)^^/^'(/' ™to the energy 
Et^, changing from (i, r, 9, (j)) to (t, r*, 9, 0) variables, and then integrating by parts 
in the variable. 

St 



(T^^)^-|-^^*^(5,V') + l^^g^ndi,^){d,^)d'f, (4.9) 
+ 2 + 

2 



a 

We can simplify the final term using integration by parts to obtain 

(dr( , ^ dr^J=]+A( dr^J=y] ^^d^* 

^^^'dV*. 

We were able to drop all boundary terms in the integration by parts, since the 
boundary terms involve bounded factors times A, which vanishes at r+, or involve 
factors of ■0 or its derivative which are assumed to decay rapidly as r — > cx) on 
surfaces of fixed t. (Merely decay such that r'-^^^tp — > or r^^^t/j — > is sufficient 
here.) Thus, we find 



By a similar argument, using the Cauchy-Schwarz inequality on (4.9 1 



ri JS' 



S2 (r2 + a2)2'^""'^' 2(r2 + a2) 



< 



r2 



< 



A 



l^'ldV 



52 (r2 + a2)3 



A2 

iV-'IdV*. 



(r2 + a2)4' 
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As a consequence of [7, equation (36)], there is the following Hardy estimate for any 
non-negative, continuous weight x which is positive at zero 



1 



(4.10) 



with the implicit constant in < depending on x, but not on r^j^ and rH,2, as long as 
< —1 and 1 < r*2- Since Q** = 1, we can apply this result to the energy for ^, 
so that 



S2 



> 



r,2 



r.i JS 



r.i JS 



-(T^^)^*^(9^^) + \gP'^{dpi,){d,i^)dh 



l^rd^M*- 



1 + 



Since A2(r2 + a^)-* < (1 + rl)-^ , 

I^TjV']*|dV< 



In JS^ 
By the same argument 

(•'"2 



r,2 



A 



< 



^ (r2 + a2)2 



□ 



Corollary 4.4 (Local decay for tp). If ip is a solution to the wave equation 'Oi'ip — 0, 
then the right-hand side of (3.321 is greater than 



Ti J-odJS 



A 1 

2 V (^^ + '^^) 



J.Z _|_ J.- 
A lr^3Af 



(19*^^1^ + iyV^I^)kA^*- (4.11) 



Proof. The right-hand side of (3.321 is 

I-T2 poo 



Ti Jr 



A2 

r.4 



1 



\dA\l + ^iv'i' + (19*^1^ + m\D ) dV 



r 



Making the substitution Tp = ■!/j(r^-|-a^)~^/^ introduces an extra factor of {r'^+a'^)~^ , 
except in the factor with the radial derivatives, where the situation is a little more 
complicated. The norm \u\2 only introduces derivatives which commute with the 
radial derivative, so we may ignore them. We will use v to denote ip with possibly 
a derivative operator acting on it. 

The sum of the radial derivative and lower-order terms has the form 

2 



A^ 



dr 



1 1 

J.2 j,2 _|_ Q,: 

A^ 

r4(r2 -t- a2) 

A2 



\vVdr 



\dri\^ -2 



;\drv\\v\ 



1 1 



(r2 -t- a2)2 (r2 + a^)^ r'^ 



liiPdr. 
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Applying Cauchy-Schwarz to the term involving vdrV and using the extra positivity 
from the original lower-order term, we have that 



A 



4 



dr 



2 



poo 1 1 



Here we have made use of the fact that + ~ . 

Making the change of variables from r to , introduces an extra factor of A/ (r^- 
a^) in the measure. Thus, the right-hand side of (3.32 (dominates 

cT2 



Ti 



These coefficients dominate the coefficients given in the statement of this theorem 
are dominated, so the desired result holds. □ 

Because we will need to work with this later to integrate by parts in the angular 
variables, we introduce the following Morawetz density with additional angular 
derivatives. Note that since the angular derivatives do not commute with §2, this 
is not quite the same as the Morawetz density evaluated on angular derivatives of 
u. 

Definition 4.5. We define the higher-order angular Morawetz densities to 



2 



(r^ + a^) r 



(r2 + a2) r4 

(Recall Oi — {©j} was defined to be the set of rotations about the coordinate axes 



in section 



24) 



4.3. K in Kerr. We now define the K vector field in analogy with the correspond- 
ing vector field in the Schwarzschild or Minkowski spacestimes. It is well known 
that this type of vector field generates a stronger energy which can be used to prove 
decay. 

Definition 4.6. We define 

K + + l)Tx + 2tr^N^dr, , 



where N is defined in (4.3 1. 

We now prove K is timelike. The essence of this argument is the same as in 
the Minkowski or Schwarzschild spacetime. The length of is very close to that 
of 9*, so (i^ -I- r^)T^ -I- 2trdr, is timelike when t^ + r1 dominates 2ir*. When 
t^ + 1"1 is comparable to 2ir* , {t^ + r^)T^ + 2tr^,dr, can become null or even slightly 
spacelike. By adding a little more T^, we have built the always timelike vector- 
field K = (i^ -I- -I- 1)T^ + 2tr^dr,. This trick of adding a little more of the 
timelike vector-field to get a globally timelike K is common in the literature for the 
Minkowski or Schwarschild spacetime. 

Lemma 4.7. There is a positive constant a, such that if \a\ < a, then K is timelike 
with respect to g. 
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Proof. Since is timelike and orthogonal to 9^. , it is sufficient to show that the 
absolute value of the norm of {t^ + r1 + 1)T^ dominates the norm of 2fr*9,-^ . 

For (t,r*) close to (0,0), the term T^^ dominates the contribution from all the 
other pieces of K, so that K is timelike in this region. We can therefore assume 
that at least one of < or |r*| is bigger than some fixed constant in the rest of this 
proof. By taking a sufficiently small, we may assume x is constant for r» bigger 
than this constant. (Alternatively, we could take K = (t^ + + C)T;^ + 2trN'^dr 
with C sufficiently large, so that this step in the argument would not impose any 
further restriction on the size of a.) 

To prove the lemma, we break the spacetime into three different regions. In the 
limit a ^ 0, the surfaces t = describe light cones. For a 7^ 0, the surfaces t = 
are no longer light cones, although give a rough approximation of the location of 
the light cones. In particular, for a small e > 0, if a is sufficiently small, the regions 
|r*| < (1 — e)t remain inside the light cone. 

Before considering different regions, wc note that 



First, we consider the region < — (1 — e)t. In this region, the difference in 

norm between and T_l is CaA^, so the absolute value of the norm of is 
bounded below by A(r^ + a^)~^(l — Ca^A), and the norm of dr, is bounded above 
by A(r2 + a^)"^! + Ca^A). The coefficient (i^ + + 1)^ dominates {2tnf + 1^. 
For a sufficiently small 



Second, we consider the region —(1 — e)t < r, < (1 — e)t. From the estimates 
on the norms of Tj_ and dt, it follows that, uniformly in the region currently under 
consideration. 




-g{dudt) =1- 



2Mr 




\{e + rl + lf9{T^,T^)\ >{e + rl + if^il ~ Ca^A) 

>((2tn)2+t2)A(l~Ca2A) 



>(2fr*)2A+ (l/2)i2A 
>(2fr*)2A(l + Ca^A) 
>{2tr,f\g{dr,,dr,)\\ 




In the region under consideration, we also have 
|<±r,| >et 



+ rl >\2tr^\ + tH'^ 
{e + rlf>{2tr^,f + eH^ 



>(l + eV8)(2tr»)2 + (eV2)i^ 



Thus, for sufficiently small a based on e 



|(i' + rl + Ifglj^, TJI >(1 + ey^){2tr,Yg{dr,,dr,) 

+ (eV2)i4|g(T„T^)|. 



Thus, K is timelike. 
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Finally, we turn to the region {l — e)t < r^. (We may also assume we're sufficiently 
far from = (0,0) that = dt.) Here, -g{T^,Ty) = -g{dt,dt) = 1 - 

2Mr/Y. > (1 - 2Afr-i)(l - Ca^r'^), and g{dr,,dr,) = AS(r2 + a'^)-^ < (1 - 
2Mr~^){\ + Ca^r^^). Following the same type of argument as we used in the first 
region, we have 

{e + rl + lf\g{T^,T^)\ 

>(1 - 2Afr-i) {{t^ + rlf{\ - Ca\-^) + + rl){\ - Ca\-^)) 

>(1 - 2Mr-'^) {{2tr^f + t^{l - Ca^)) 

>(1 - 2Mr-^){2tr^Y{l + CaV^^) 

>{2tr.^fg{dr,,dr,). 

□ 

We now show that the K momentum dominates the momentum and that, on 
hypersurfaces of constant t, the energy assciated with (Jk is small relative to the K 
energy, so that we can freely add or drop the energy associated with q^^ when we 
find it convenient to do so. 

Lemma 4.8 (K energy). There are positive constants a and Cke, such that if 
\a\ < a, and if is smooth, then : 

(1) (estimate on the energy) 

CkeVI, >\{t + r,)(T^ + drj^\^ + \{t- r,)(T^ - a,J^|' 

+ (t' + rl + i)(/i"''(a„^)(9^^) + v^p^). 

As a consequence, for |r»| < (3/4)t, 

CKE-Pk > (l + i')^T,- 

(2) The contribution to the energy from gK is small in the sense that 
\EgAt)\ = \Q'^qi^{dp^)^ - \Q"'{^|3ql^^''\d'^i., <a^CKEEK{t). 

Proof. This follows the same ideas as the previous lemma. We start by computing 
the energy, 

=^N-\T^'tJj){t^ + + l){Tx^P) + {T±i>)2tr,N^{dr,'ijj) 

-^-(t^ + rl + l){N-\T^^jf ~ {drM - h'^\d^^){dfii,) - Vi^^) 
= \if + rl + \){N-\T:^i,f + {drM + h''l'{d^i,){dpi^) + V^^) (4.12a) 

+ {t" + rl + 1)N-\T^^){u;hX ~ ^iMd^i^) (4.12b) 

+ 2tr.,N'^{T:i_ij){dr,^). (4.12c) 

We now seek to use ( |4.12a ) to control the remaining terms. 

First we consider the region — t| < t/4. In this region, |tr*(A^^ ~ 1)1 ^ 
and \{t^ +rl + \){N'^ - 1)| < a, so that the term (Tj^t/))^ + (9^.-0)^ dominates any 
terms arising from replacing all the factors of TV by 1. Similarly, since ujhX ~ "^J- 
a/S./{r'^ + a^Y , the term (4. 12b I is dominated by a small multiple of (4. 12a I. In 
fact, because of the weights, the lohX — aA/{r^ + a^Y term can be controlled 

by a 

+ rl + ^)(iV-2(T^V^)2 + [drM + h^^{d^^){dpi^) + V^^) 
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Thus, in this region, 



t^ + rl + 



> 



|-2in(Tj_7/.)(5^.V') 

+ (f + rl + l)(h''^{do.i;){dpi,) + Vi>^). 

In the region ||r*| — t\ > t/4, the dominant term [t^ + r^)((Tj^?/')^ + idr,ip)'^) 
domin ates (t^ + rl){N-^ - 1)((T_lV')^ + (dr^'ipf) and the terms in ( |4.12b| and 
(4.12c I by a fixed factor greater than 1, so that the desired estimate holds. The 
control over follows from this. 

We now turn to the terms arising from gK- These terms are 



<t\N^ ~l\\T^i;\\4,\ + ^\N^ 



(r2 + a2)' 



r^2 + a2(T^^)2 



(r2 + a2)' 



The first term is controlled by "P-^, and the second, by . The decay of the final 
term is too slow to be directly estimated by the energy, but after integrating along 
a hypersurface of constant t and applying the Hardy estimate (4.101, the third term 
can be estimated by Et . Thus, 



1 X 



St 



□ 



Lemma 4.9 (Bulk term for K). If ip is a solution to the wave equation (4.11 and 
tp = + a^^, then 



-d. 



gK),3 



< 



{i + ty^v^ 



for r, < -t/2 



and, in the region — </2 < < t/2 
1 



((1 + log(2 + 3 for r, > t/2 

t/2, 



<t{5VAm 

+ t{5VA)P\,^P] 
+ \a\t\5VAm 

K,3- 



\a\t-'Vi 



Proof. We begin by considering 

K =(t2 + rl + 
Tk =dt + uj-i^d^ 



l)Tvi + 2tr^N'^dr,, 



and derive conditions on wk- Near the end of this proof we will conclude that 
Tk = is permissible. For simiplicity, for most of this proof, we use v to denote 
a second-order symmetry operator acting on ip, and consider the deformation of 
P(K,gK)[^]- Similarly, unless otherwise specified, and T-'k are understood to be 
evaluated on v. 
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The bulk term we want to control is 

(4.13a) 

-qKg'"^{dc.i){d0v) (4.13b) 
+ (^l{d^K'')V + iK^d^V) - qkV^ v\ (4.13c) 

+ \{d0g''^d^q^)v\ (4.13d) 
With an eye to computing the derivative terms, we first compute 
d-^KT' =2t + 2iiV2 + 2tr^{dr,N^) 

=4tN^ + 2t{l - + r.dr.N'^). 
The terms involving the derivatives and arising from the vector field are 

- {d,v)2r4Tj^v) - {d,v){t' + rl + l){dr,wj^){d4,v) 

- {d,v)2tN^{d.v) - {d,vf2tr.{dr,N^) 

- h"'^{dcv){t'' + rl + l){d0U>K){d^v) - h''^2tr4d0N^){dr,v) 

- (2tiV2 + t{l - + r,dr,N'^))N-'^{T^df 
+ {2tN^ + t{l -N^ + r,dr,N^)){d,vf 
+ {2tN^ + t{l -N^ + r,dr,N^))h"l^{d^v){d0v) 

+ tu{dr,h'^^){do,v){d0V). 

Expanding Tk as Tj_ + (wk — u)±)d^ and grouping terms, we find 

{d„v){d0v) (^-g'^^d^^" + ^(a^K'^)^"'' + ]^vi'^{d^g"^'] 

= - {T^i)H{2{l - iV-2) + iV-2(l -N^ + r.dr,N^) + r.dr,N-^) 

(4.14a) 

+ {T^v){d^v){2t)N-^{ujj<: - u^) (4.14b) 

- {d,v)H{2n{dr,N^) -{1-N^ + ndr,N^)) (4.14c) 

- {d,v){d^v){{t'' + rl + 1)(9^.wk) + 2n(a;K - (4.14d) 

- {d,v){d0v)h''^2tnid0N^) (4.14e) 

- idc,v)id0v)t{-2N'^h"'^ - nOr^h"^) -{1-N'^ + r^dr,N^)h°'^) 

(4.14f) 

- {d^v){d4,v)h"^{t^ + + l){dpcvK). (4.14g) 

In considering these calculations, it is important to note that the worst terms, 
involving {T±v){d^v) are exactly cancelled. To force this cancellation, we have 
included TVMn K = (i^ + + 1)T;^ + 2N^trJ)r, instead of using [t^ + rl + l)dt + 
2tr^dr,, which is a more straight-forward analogue of the vector field used in the 
Minkowski and Schwarzschild spacetimes. We now analyse the remaining terms. 
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For the terms (4. 14a I , ( 4. 14c I , and ( 4. 14f I , which we consider to be hke derivatives 



squared times t times coefficients, we look at just the coefficients. We also include 
the contribution from qnG"^ ■ The coefficients of {T±ip)'^ and {dr,ip)'^ are 

Coefficient from ^Ula^ ^ - N-'^{N^ - 1 + r.,dr,N'^) - N^r^dr^N^^ + {qK/t)N-^ 



=0. 



H^'-i-(gKA)) 



Coefficient from ( |4l4c| = - 2r.,dr,N^ + (1 - + r^dr,N^) 

= - + l~r,dr,N'^ - (qK/t) 

= - 2{N^ - 1) + N^r^dr,N^^ 

=a^ sin^ eO{r^/^, r^^ log r). 

Thus, for < ~t/2, with the t{dr,v)'^ factor, the term is of the form a^fr*A(9, „u)^ < 
a'^t-'^{dr,v)'^ < t-'^V^ . For -t/2 < < t/2, the term is bounded by a^r 



For < t/2, the term is bounded by a'^{t-^ log t)V^ < a^{t-'^ log t)Vl^. 
The angular derivatives are given by 

+ (2(7V2 - 1) + (1 - N^) - u{dr,N^) ~ {q^/t))h''^. 



Coefficient from (4.14fl 



the first two terms of the coefficient of (4.14fl are of a familiar form from the 



Schwarzschild spacetime. Outside a compact set near the orbiting nuU-geodesics, 
they have a good sign, i.e. are negative, and decay like A near the event horizon 
and not slower than r~'^logr as r ^ oo. The remaining terms involve only factors 
involving either {N'^ ~ 1) or r^,dr,N^'^ multiplied by the potential Vq in h. Thus, all 
the terms should decay at least like log(A)A^ near the horizon and at least as fast 
as r~'^ as r ^ oo. In addition, these all have a factor of a on them. Thus, outside 
a slightly larger compact set, the negativity of the first two terms is sufficient to 
control the rest. Inside that compact set, the sum, with the t{da'4^){dpip) factor 
included, can be controlled by 

We now turn to showing the remaining terms are small. The terms we wish to 
control all involve mixed derivatives. We first estimate these away from photon 
orbits. We break the remainder of the spacetime into the three regions /, //a, and 
///, where < —t/2, where |r*| < t/2 but |r*| > C to avoid the photon sphere, 
and where jr*! > t/2 respectively. Region Ilh will refer to the region |r*| < C. 
To estimate the remaining terms, we will apply the Cauchy-Schwarz inequality and 
then use the following estimates in the regions /, Ila, and ///, 



1 



1 



(r2 + a2) 



respectively. Typically, in region /, we will seek an additional A*^ decay, whcih we 
can estimate by i^^. Similarly, in region ///, or r decay can be traded for t 
decay. Since we are applying the Cauchy-Schwarz inequality, these estimates limit 
the asymptotic behaviour of the coefficients of the mixed derivatives. This imposes 
restrictions on the behaviour of cjk and its derivatives. Thus, in these calcula- 
tions, we find ourselves having certain coefficients, needing them to satisfy certain 
asymptotics so that we can apply the relevant estimates, and then requiring certain 
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conditions on smd its derivatives. We call this "need-have-require" analysis. 
When there is no condition which we require, we simply say the required condition 
is "nothing" . 



The first term, coming from (4.14b I, is 

In the three regions of interest, the asymptotics we have, those we need, and the 
conditions we require are 



Region 


Have 


Need 


Require 


/ 

Ila 
III 


t{ujK - ^^±) 
t{uJvL - UJA_) 


|a|iO(A,r-3) 
r-2 


^K-^i = |a|0(A,r"3) 



The second term, from (4.14d), is bounded by the sum of two terms 
\{d,v){d^v){{e + rl + l){dr,L0-K) + 2n(cjK - 0J^))\ 



< 



\d,i\\d^i\ {{f + rl + 1)|9..wk| + 2|n(u;K - 



We perform a have-need-require analysis on each of these. For the first of these 
terms, in region Ila, we will not estimate the term by t{5VA)['4'\ but by af^((57'^ )[?/;], 
which means we "have" an extra factor of t. 



Region 


Have 


Need 


Require 


/ 

Ila 
III 


if + rl + l)dr,uj-K 
{f + rl + l)dr,uj-K 

{t' + rl + l)dr,UJ-K 


Al/2+e 

ai2o(A,r-5/2) 


\dr,LO-K\ < Al/2+^ 

i9r.WK = |a|0(A,r-5/2) 


For the second, we again use the standard estimates. 


Region 


Have 


Need 


Require 


/ 

Ila 
III 


r^^K - 


Al/2+e 

|a|iO(A,r-5/2) 


c^K-c^± = |a|0(A,r-5/2) 



The next term from ( 4.14e[ ) is estimate similarly. These are 
\id,v)idf3v)h'"'2tr4dpN^)\ <\d,v\mV§t\r,\ 

and we require 



Region 


Have 


Need 


Require 


/ 






nothing 


Ila 




iO(A,r-5/2) 


nothing 


III 






nothing 



For the final term, from (4.14gl, we again need to use at'^{SVA)[ip] to dominate 



Region 


Have 


Need 


Require 


/ 








Al+e 


ly^Ki < A^ 


Ila 






lyt^Ki 


ai2o(A,r-3) 




III 






IV^kI 
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We turn to controlling these terms in region lib, near r — 3M. If we assume 
that is constant in a neighbourhood of the photon orbits, then the worst terms 



involving t , coming from (4.14dl and (4.14gl, vanish in this region. We will further 



assume that in this region t^K — uj± is bounded by a constant times a, which bounds 



the term coming from (4. 14b I. Applying Cauchy-Schwartz, all the other terms are 
bounded by 



tlr 



atlr 



z3M 



\dtv\ 



By taking the most restrictive conditions, we are left with the following require- 
ments: 

uJK-oj± =|a|0(A,r-3), 

dr,cuK=\a\0{A,r-^), 
yc.K=|a|0(A^r'3)^ 

and that wk is constant in a neighbourhood of the photon orbits, and that cjk — ^± 
is bounded by a constant times a in that neighbourhood. 

There is a wide range of possibilities for choice of wk- In particular, our choice 

satisfies all the necessary conditions. 

The remaining terms to be estimated, which come from (4.13c )-(4. 13d), are given 
by times the coefficient 

^{d^K-')V + ^K^d^V - qkV + ^{daG^^dpqK) = + tN\2V + r,dr,V) 

+ t{N^ -l + r,dr,N)V 

- (N' - l)qK ~ lidlqK) 



As with the angular derivative terms in (4.14fl, the first term is familiar from the 



Scwarzschild case and gives a weight which is positive in a compact region, but 
which is negative, and hence trivially bounded, outside this region. The behaviour 
outside this region is — tA log |A| approaching the event horizon and —tr~^ log r 
approaching the horizon. The other terms are all a^O{A,r~'^) at worst, so they are 
easily dominated by the first term outside a slightly larger compact region. Thus, 
the contribution involving without derivatves is easily controlled by t{SVA)[4']- 

□ 



4.4. Stronger, light-cone localised Morawetz estimates. In this section, we 
prove a stronger Morawetz estimate to gain factors of t inside the light cone, |r* | < t. 
We do this by using a Morawetz vector field with additional factors of t and localising 
inside the light cone. 

Definition 4.10. Let xlc : M ^ M &e a smooth, even Junction which is identically 
1 for \x\ < 1/2, identically fori \x\ > 3/4, which is weakly decreasing for x > 0, 



The parameters 1/2 and 3/4 are not important. The important thing is the ordering < 
1/2 < 3/4 < 1. 
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and which is weakly increasing for x < 0. If Xlc written without an argument, 
it understood to mean 



Xlc =Xlc 



(t) 



The stronger, light-cone localised Morawetz vector field (with strength p) is defined 
to he 



Af =tPxLcA'^, 

-ab ,ri ab 



Note that 



—A— Xlc, 

-ab ab 

1a,o =<1aXlc- 

Lemma 4.11. There are positive constant a and C, such that, if \a\ < a, if ip is a 
solution to the wave equation ~ 0. and if = {r^ + a'^)^^'^^, then for t2 > ti > 
and p e 

r I tPxLc{5VAm<i'i^<{^ + hY~^E^M{t2) + {i + hY-^E^M{h) 

''\l + tY-^E^^S]{t)dt. 

Proof. We begin by relating the divergence of the momentum of the stronger hght- 
cone locahsed Morawetz vector field to that of the original Morawetz vector-field. 
We denote the right-hand side of (3.6 1 by W [■(/;]. The divergence is 

-^^{^iVA,)^-^o.{^J^tPxl.cVA) 

XLC -dc{f^V%) + V^id^XLctn, 

- r^ida.tPxLc) - tP XLC -dciKBAj + Ba-ji))- 
We now move the localisation onto the solution itself to find 
i^W[xLcV'] 

- V2{datPXLc) - tP XLC- do.WA;l + Ba-Ji)) 

- ^"^TTT-^ Y.^2x\!^^d^rSJ^){drx\!c)S^i^ + {drx\!c)\SJ^^)), 

^^^d„(f,tPxLc{V2 + B%, + B%u)) 
-iV2 + B%, + B%uKdo.tPxLc) 
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We now estimate the last term, using X2 as a smooth locaUsation which is 1 on the 
support of xlc but still supported well inside the light cone, and using <p to denote 



»f A2 
<t^X2 (l^{dr4>f + 



Integrating this over a hypersurface of constant t, 



I i^X2^T^^((a.^)(9.XLc)V'+ (5.XLc)V')dV 



Since V\, V\, B%j, B^.j, B^.n, and B^.n are all bounded by ^^-^[§2^] + |S2V'P, 
which is bounded by [S2V']; at least in an integrated sense over a portion of the 
hypersurface inside the light cone, it follows that 



+ tp- 



/■(3/4)t r 
J-{3/4)t JS^ 



Now applying an argument similar to the one earlier in this proof, it follows that 



/ XLct^WlV-jdV < / -5„(Mt^XLc(?'3 + SS;i + B2.„))dV 

J-(3/4)t JS^ 
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To finish, we integrate in time, and then apply the Hardy estimate as in lemma [3.12| 
and estimate the lower-order norms as in lemma [3.13[ This gives 

rt2 /•(3/4)t 



Jti J-(3/4)iJS2 
ft2 ri3/4)t f 

+ / t"-' / / njS2V^]dV* 

Jti J -{3/4)1 J 

<ir'^K[§2V'](i2) + tr^EK[S2mi) 



□ 

4.5. Closing the K estimate. In this section, we use the stronger, light-cone 
localised Morawetz estimate to control the growth of the K energy. 

Theorem 4.12. There are positive constants a, C , and C such that, if \a\ < a 
and if ip is a solution of the transformed wave equation (4.1 1, then Vt > 



EK,3it) <Cil + 0^''"' (^K,3(0) + ^T„5(0) 



Proof. Essentially, the stronger Morawetz estimate, lemma [4.11| controls the growth 
of the K energy, and the K energy controls the boundary terms in the stronger 
Morawetz estimate. Thus, a standard Gronwall's argument closes the estimate. In 
this proof, an integral without limits typically refers to integration over (r^,9,(j)) € 
R X 5*2. 

First, we can bound the growth of the K energy, using lemma |4?9| by 

rt2 



^(K,,K).3(i2) - ^(K,,K),3(tl) < / ' ^T„3(t)(l + T^dt (4.15) 



tl 

4 



t{5VA)P\Md''ti* 

e{5VAmd''^Ji, 

+ \a\ [ {l + t)-'EK,3dt 

+ / \l + t)-^logi2 + t)EK,3dt. 
Jti 

Since Et^^s < E-k,3 and (1 + t)-"^ < (1 + i)"2log(2 + t), the first term can be 
absorbed into the last. 
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Before applying the Morawetz estimate, we apply integration by parts in the 



angular derivatives to the term in (4.15) and apply the Cauchy-Schwarz inequality 
to find 

3/4, ^ I, \ 1/4 



Now applying the integration in the remaining variables and the Cauchy-Schwarz 
inequality one more time, we find 



t2 

Iti Js 

From this, we can bound the growth of i?K,3 by 

rt2 



((57'A)[Ot,^]dV* 



^K,3(t2) - ^K,3(ii) <r I i^/'(<57'A)[^]dV* 

(^7'A)[Ot,^]dV 



i«i r I i'('5^A)[^]dv* 

Jti JT.t 

\a\ I (l + i)-i^K,3di 
/ V + i)"'log(2 + t)^K.3di. 



In the first integral, with the aim of applying lemma [4.1H we use the convergence 
of t^/^~^dt to get a bounded piece plus what we hope to show is a small piece. 
To do this, we choose a time T such that /^(l + t)~^/^dt is small relative to the 
absolute constants implicit in <. The first integral we break into an integral from 
ti to T and then from T to is- In the first part, we replace (1 + 1)^/'^ by (1 + T)'*/^, 
which is also an absolute constant. Thus, the first integral is bounded by 

*2 !■ poo !■ 

/ t^/^5rAmd^^i*<E:,.T^+ / t''^5VAmd^^i*. 

ti -/St JT J^t 

Similarly, in the last integral, we want to use the integrability of J t^^logtdt 
and the bound En (1 + ^)^^Tj^, so that we can take T sufficiently large that 
log tdt is small, and obtain the bound 

/ ' ^K,3(l + log(2 + t)dt <Et^.3 + f ' ^K,3(l + t)-^ log(2 + t)dt. 

Jti JT 
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We can also apply the original Morawetz estimate, lemma |3.13| with an extra 
four derivatives, to control the integral of ((57'^)[0^, -i/;]- This leaves us. 



Jt 

+ |a| r t\5VAmd''^Ji, 



2 1 



£^K,3(l + t)"'log(2 + t)dt. 



Using the stronger Morawetz lemma 4.11 we can bound the integral of the 
Morawetz terms by, 



Ej^Ah) " Ej^Ah) <(1 + t2)"'/'^K,3(i2) + (1 + r)-2/3^K,3m 
+ JJ{l + t)-^/''ET^,3mt 

+ WlEnAh) + WlEnAti) 

+ \a\ f\l + t)-'Eji,3{t)<it 
Jti 

EK,3{^ + t)-^log{2 + t)dt., 



(4.16a) 

(4.16b) 

(4.16c) 
(4.16d) 

(4.16e) 
(4.16f) 



Let 



f{t) = sup Ek,3{t), 

rG[ti,t] 

Kq —EK^3{ti) + Et j{ti). 



We can replace E'k.s on the left of (4. 16a I by /. We can make a similar substi- 

1 1 1 1 1 1 1 -1 2/3 

tution in (4.16al, ( 4.16b[ ), (4.16d), and (4.16f), and then use the smallness of ^2 i 
gj^^ ^"^logt to absorb all those terms into the left-hand 
side, with only a small loss. Note that by the bounded energy result, one has that 
Et ,7{t) < Kq uniformly in t. Thus, we are left with 



f{t2)<Ko + \a\ / {l + t)-'f{t)dt. 
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We now write this as an integral inequality to apply the standard Gronwall's in- 
equality techniques 

f{t)-C\a\ f (l + r)-V(T)dT<Xo, 



il + t)F'{t)~C\a\F{t) <Ko, 
d_ 
dt 



Fit) <^Xo(l + t)^l'^l 
\a\ 

fit) <Xo(l + t)'^l''l. 
In terms of the original quantities of interest, this gives 



□ 

We also want to evaluate the K energy on surfaces other than those of constant 
t. The following lemma allows us to control the K energy on such surfaces. Al- 
though this lemma is valid for any surface, it is only interesting on time- like or null 
hypersurfaces on which the energy should be positive. 

Theorem 4.13 (K bounds for other surfaces). Consider the almost-null cones 
given by 

Cr = {t~\r,\ = T}. 

There are positive constants a, C , and C , such that, ij\a\ < a, and if'tp is a solution 



of the transformed wave equation (4.1 1, then for any hypersurface S, which need not 
be in the futute ofT.Q, if T. can be represented as a graph over (r*,0, ^) g M x S*^ 
which lies beneath Cr and above —Cr, then 

^(k,,k).3(S) <Cmax{l,r}^>l (£;k,3(0) + £;t„7(0)) 
Proof. In this proof we will exploit the time sy mme try of the previous results. For 



simplicity let Kq — £'k.3(0) + £'t;^,7(0). Lemma 4.9 allows us to estimate the change 
in -©(k.ijk) from {t = 0} to S by the four-dimensional integral of da{^jLP"^ ^k))' 
This integral can be broken into three pieces: where r^ < — |t|/2, where — |t|/2 < 
r* < |i|/2, and where \t\/2 < r*. 

We first consider the middle region, where — 1<|/2 < < \t\/2. The four- 
dimensional volume inside this region and beneath S is a subset of the volume 



inside this region where \t\ < 2t. Thus, from the proof of theorem 4.12 the integral 
of the absolute value of ^^^da{^iP'^j^ ^^■^) is bounded by C(l -I- t)^°^q. 

We now consider the other two regions. First, consider the fax region, \t\/2 < r^. 
In this region, |M"'5„(/iPfK,?K))l bounded by log(2 + |t|)(l + If we 

now foliate the region beneath S by hypersurfaces of constant t, we can bound the 
four-dimensional integral of fi^^da{fJ,V"j^ ^^-j) by the integral in time of the integral 
over each of the leaves of the foliation. On each leaf, the integral is bounded by 
log(2 + |i|)(l + |i|)-2^(K,5K)W < log(2 + \t\){l + \t\)^\-\-^Ko. Thus, the four- 
dimensional integral is bounded by Kq f log(2-|- \t\){l + \t\)'^^°-^^^dt < Kq. The near 
region, < — |t|/2 can be handled similarly. 

Thus, the four-dimensional integral is bounded by Ct'~^^°'^Kq. Hence, so is 

^^(K,gK)(^)- n 
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5. POINTWISE DECAY ESTIMATES 

5.1. Decay in stationary regions. In this subsection, we prove decay in regions 
of fixed r bounded away from the event horizon, < ri < r < r2. Since these 
are preserved by the flow of the stationary, KiUing field, dt, we refer to these as 
stationary regions. 

The essence of the proof of this decay result is that the bound on i5K,3 gives a 
bound on times the square of a local norm, which controls the square of the 
field, through a Sobolev estimate. Thus, |V'| decays Hke t^^ with a small loss 
from the growth of the K energy. 

Theorem 5.1. There are positive constants a and C such that, given r_|_ < ri < r2, 
there is a constant, Cri.r2; such that, if \a\ < a, and if ip is a solution to the wave 
equation dip = 0, then > 0, r e [ri,r2], {9,4)) € S'^, there is the estimate 

mt,r,e,<t>)\ < a,.r.max{l,i}-'+^'l"l (^k,3(0)^/^ + i?T, jW^/^) . 

Proof. In this proof, a constant of the form C^^ refers to a constant which depends 
only on ri and r2. As is common in analysis, the value of this constant may vary 
from line to line. Since ?'+ < ri < r2 < oo, there are corresponding values 
— oo < r*i < r^,2<^oo■ 

for fixed (i, r,) and all {9, cj)) £ 



From lemma 



2.2 



m,r,,9,4>)\'< / IV-bdV 

Given any smooth w : M — > M, by a one-dimensional Sobolev estimate (or simply 
through the fundamental theorem of calculus, localisation, and Cauchy-Schwarz), 

for r.^,l < <r^^,2 

Kn)P<C;,,, r^\drM^ + \v\^dr,. 

Applying this result with v being a symmetry operator applied to "0, then summing 
over the symmetry operators, and then integrating over S"^ , we find 

IV^bdV < C;,r. / n \drM + IV^IsdV*. 

The left-hand side dominates Since the range of r is fixed, there is a lower 

bound on y, so that 

r r^*2 

|2 ^ I I \ ^/.|2 I t/L/,|2j3, 



S2 



For sufficiently large t, the region [r'*i,r*2] x is inside |r*| < t/2, so that the 
right-hand side can be bounded by the K energy, and 

m,r,,9,4)\^ <C;^,r,t-'^K.3(t) 



<c;^,,,t-2-^'^ (^K,3(o) + ^T,j(o) 



Since (r^ -l-a^)^/^ is uniformly bounded above and below on (ri, r2) we may replace 
tp{t, r^,9, (j)) by 'ip{t, r^,9, cj)) to get the desired result. 

For small t, the result holds from a similar argument using the energy bound. □ 

5.2. Near Decay. To study the wave equation near the event horizon, it is conve- 
nient to use u+ and u_, which were defined in the introduction, and to use 

The functions {u^,r,9,(f>^,) are known to form a coordinate system which can be 
extended to an open set containing the future event horizon [57] . 
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Theorem 5.2 (Near Decay). There are positive constants a, C and C , such that 



if \a\ < a and if tp is a solution of the wave equation (4.1 1, then for all t > r*, 
r < MI, (9, (/)) e S^, 

|V(t,r,0,0)| <Cmax{l,u+}-i+^''^ (^k,5(0)^/' + i?T„9(0)i/2 + 4.,^ ,3(0)^/') . 



Proof. Near the horizon, it is convenient to rewrite the wave equation, (4.1), for ip 
as 

+ A{Q + dl)4' 

2 \ _ 



- 

= -{dt + ujHd^y^p + + - V^, 

where Z can be expanded as Z = Z-Sa and the Z- vanish hnearly in r — r+. 
Equivalently, the functions Z-/ A are bounded for r G [r+ , 2>M] . 

To prove decay, it is convenient to use u_, 61, as a coordinate system in 
the exterior region. The coordinate derivatives in this coordinate system can be 
related to the earher (i,r*,0, 0) coordinate derivatives. The relations that we are 
interested in are 

d± = 9u± =^ {dt + i^ud^ ± 9r. J , 
de =de, 

The relation dg — dg is understood to mean that the ^-coordinate derivatives in 
the two systems generate the same vector field. We can, without ambiguity, write 
formulas involving coordinate derivatives from both systems. In particular, we can 



write the wave equation, (4.11, as 



= - 49+a_?A + Z^P-V'tp 

with Z defined in terms of dt, dg, and c?^ as above. 
To prove decay, we use the one-form 



^ — i^d-tpj du-. 

For a given {6, cj)^), let 

n^{{u'^,u'_,0,(j),)\u'_ e e [-■"'-,"+]} 

={{t'r[,9',(l)')\t' >0,t + r,<t'-ri <t~r,,t' -r', <t + r,}. 

It is convenient to first consider u+ > 1. We apply Stokes' theorem to £, on 
and in m_ 0, 0^) coordinates to get 



dC = / 

n Jon 

(9+9_'0)du+du_ =ip{t, r^,6,(j)~- iVHt) — + r^,0,9,(j) — uJnit + r*)) 

{d^i>)dr,, (5.1) 

c 

where C is the curve with t' = 0, e [— u_, and {6', 0') — {9, (f)). 



n 
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Figure 1. The region il in the {u'^,u'_) plane with vertices at 
{—U-,U-), (— and «_)_). In the (f, r*) plane, 
these vertives correspond to {t,r^), {0,—t + r*), {0,t + r*), and 

(t + n,o). 



The function evaluated at zt+, 9, 0,) corresponds to the value at (t', (j)') 
0, 9, — UJHU+), so it can be estimated by 

The integral over C corresponds to an integral in the hypersurface t = with 
e [~U-,—u+] C (— oo, —u+]. Since 9_ is in the span of Tj_, dr^, and Ad^ with 
smooth and uniformally bounded coefficients when r G [r_|_,3M], it follows that 

|(A-V2a_)^|2 =A-V_Vi|2 < A-ip*,^. 

Since Tj_ is parallel to the normal, ns^ , it follows that 

|(A-V2a_)Vi|^<A-V2p*^^. 

Therefore, 



The second integral can be bounded by 



/ ^ |(A-i/25_)V;pAV2dr, < [ f * |(A-i/25_)Vi|^Ai/2d 



< / / ,3dV* 



S2 J-OO 



<i?n.„3(0). 

Since A decays exponentially in . which is faster than any polynomial 
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Finally, to estimate the integral over CI, we break the integral into two regions 

A=nn{n > -t/2}, 
B=nn{n < -t/2}. 

We first prove a preliminary estimate in region A. For —t/2 < r* < 0, since 
A/(r^ + a^) is increasing, it follows that, in (t,r^,9,(j)) coordinates 

Integrating along a curve of constant {t, 6, 4>), we have 
A ~ 



+ t2V'(i, 0,^,0)2 

f° 1 ~ - 

<^K,3(t) 

<f\'^\ (^K,3(0) + ^T„7(0)). 

Therefore, 

Returning to the argument involving Stokes' theorem, the integral over CI can be 
estimated by 

/ (9+9_'0)du+du_ = / (^V' ~ l^V')d'U+dw_ 
in in 

=2 / (^V' - VVi)dMr* 
in 

<2sup|AV^Vi|2 / A^/^dMr^ 

+ 2 (J^ AdMr*^ (/si (^•^'^^' ^ ^^'^^') ■ 

We now control each of these terms. We start with the first. By the prelimi- 
nary estimate in region A, there is the estimate A'^/'*|'(/'| ^ i^"'^^'^'"' (-E'k,3(0)"'^/^ + 
£'Tx,7(0)^^^). After applying up to two more derivatives, and noting i' ~ in A, 
we have that the supremum term decays like 

sup |aV4^|2 < u;^+^i«I(£;k,5(0)^/^ + i;T„9(o)^/^). 

A 

Since A d CI has \il — 2u+| < \r'^\, the integral in /^A'^/^dMr* is bounded by 
2A^/^r*dr*, which is bounded, since A decays exponentially. 
We now turn to the integrals over B. Since 2, and V decay linearly in r — 
the integral / A~^((2^V)^ + (VV')^)dMr* is bounded by a multiple of A times the 
square of up to two derivatives of V^. Therefore, after applying two more angular 
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derivatives and a spherical Sobolev estimate to restrict to B instead of a four- 
dimensional region, this integral can be controlled by the Morawez estimate 

^ i (^{Zi,)^ + dMn < i?T,,5(0). 

Finally, the integral /g Adtdr* is bounded by an integral of the form J Ar^,dr.^,, but 
now the upper bound on is ~ since we are restricted to the region B. Thus, 
Jg Adtdr^ decays faster than any polynomial in . Combining the estimates in 
regions A and B, we have 



Thus, from (5.1 1, the solution -0 decays like 



This gives a bound on at (t, r^,,6,(j) — LOnt). Since the same decay rate holds at 
all points on the sphere, the desired result holds. 

For u+ e [0,1], the same argument holds to give the boundedness of ip. For 
u+ < 0, a similar argument holds, except that f2 is a triangle in the (u_|_,m_) plane, 
with both endpoints on the surface t = 0. The boundedness of the initial data on 
t = in the region r < 3M, then gives boundedness for < 0. 

Since ^l> and 'ijj are uniformly equivalent for r < 3M, this gives the desired result. 

□ 

5.3. Far decay. We now prove decay in the far region, for r > r^, with particular 
attention to the behaviour as r ^ t and r — > oo. 

To study the far region, we introduce the radial coordinate, y, and almost null 
coordinates, v±, 

y= hdr[, 



(r2 + a2)2 



v± =t ± y. 

Since, for large r*, /i ~ 1 — a^Cr^^, the coordinate y differs from by at most a 
constant. Similarly v± differs from u± — t ± r^, by at most a fixed constant. By 
direct computation, one finds that the length of the one- form 

dv± =dt ± hdr^, 

is given by 

-—- — ^g{dv±,dv±) ^g°''^{dv±)a{dv±)f3, 

g'^l'{dv±)o.{dv±)p ^a'ism^e- 2)Vq < -g^Vq < 0. 

Thus, the surfaces of constant u_ (or v+) are timelike. Although we have chosen a 
particularly convenient form for the factor h, the only property that we use is that 
/i ~ 1 — Cr~2, so that the surfaces are timelike, with a fixed rate at which they 
degenerate with respect to the T_l and dr, basis, and that the constant in this rate 
is sufficiently large that < C. We denote the surfaces of constant V- by 
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Unless otherwise stated, we will restrict these to the region where t > niax(2y, 0). 
Outside this region, we extend the surface as a surface of constant t. It is also 
convenient to introduce the almost null vectors 



^± =Tj_±h-^dr.^. 



These satisfy = = ^^'^ — 1 = C-"*^-- Before proving decay 

in the far region, we control the K energy on almost null hypersurfaces on which 
t — = 0(1) and r — > oo. In the standard terminology, which we explain in the 
appendix, these almost null surfaces would be described as crossing null infinity. 

Lemma 5.3. There are positive constants a, C , and C such that, if \a\ < a, and 
if ip is sufficiently smooth, then 

cE^iKJ > f K + 1) f + dv+d'^. 



Furthermore, if ^ is a solution to the wave equation, (4.11, then 

EndKj <Cv^'^''^ (^K,3(So) + ^T„7(So) 

Proof. It is convenient to introduce the pair of null vectors 
The wave equation can be written as 

For any vector in the span of Tj^ and dr, , 
the vector can be written as 

X± =1 ("iVX-L ±X' 



2 

The vectors 

B± =7V"^T_L T hdr, 
have null components (ie, components with respect to the and LJ) 



We're particularly interested in the components of _B_, which can be estimated by 

BZ >a^VQ > a^r^^. 
Similarly, there is a decomposition of K as 

K =K+L+ + K L - {t^ + rl + \)oj^d^, 

((^^ + + 1)^ ± ^tr^N' 
= ^{t±r,Y + ^±tr,N{N -I). 
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For r, > tjl and a small 



|(t2 + r2 + l)c.^|<^, 
r 

K± >t;| + 1. 



The functions {v-,r^,0,(j)) can be used as a coordinate system. (In fact, the 
coordinates {v- ,r,6,(j)^) can be used as a coordinate system extending through the 
(future) event horizon.) Since the new coordinates satisfy 

dv- A dr* A fjA6 A d^ = dt A dr* A fjA6 A d^, 

the energy associated with any vector-field X is the integral over the 3-surfacc of 
the contraction of the associated momentum (vector) against the 4-dimensional 
measure 

/ 'Px(df A dn A /id6' A d(/)) = / (dt;- A dr, A yud6l A d^) 

>/e+_ J^t_ 

= [ ^^-dV- 

In {v- ,rt.,6, 4>) coordinates, the energy associated with K on a surface of constant 
V- is the integral of 

+ i<5f (dT;_)0K" (^-iL+4,){L^^P) + h'^(\dj){dp,p)) 

+ ((B^)(K-) + (SZ)(K+) - ^(dt;_)„K") (L+^)(L_Vi) 

+ (BZ)(K-)(L_Vi)^ 
+ ^(dt;_)„K« (h"f'{^„i,){^0^^>) + Vi>^) 
+ \a\Oir-^)\B^d„i,\\d4,^\. 
By direct computation, 

(dt;_)„K« =2 {{Bl){K-) + (BZ)(K+)) , 

so that 



^2 



Here, the term |a|O(r~^)|i?_^/'||90'!/'| has been absorbed into the remaining terms 
with a small loss in the constants. 
Since 

L± =N-^dt + Nw^d^±dr,, 

it follows that 
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The difference between ^_|_ and can be estimated by 

^+ - L+ =(1 - N-^)dt + (1 - h-^)dr, + N-^uJ±d^, 



A similar estimate holds for ^_ and L_. For > t/2, the desired derivative can be 
controlled by 



+ {vl + 1) (-2(L+Vi)((C+ - + ((€+ - L+)i^f) 

<ivl + l){L+i>)^ 

+ {vl + l) ((e+-L+)V')' 
<K + 1)(M)' 



Thus, the energy on a surface of constant V-, associated with the vector K is 



At this point in the argument, we want to show that the energy associated with 
gK is negligible so that the slow growth of the (K,(7k) energy implies that the K 
energy grows (at most) slowly. The energy associated with the scalar gx is 



If E+_ is understood to be restricted to the region r* > t/2, then 
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The integrand in the first integral is controlled by , so that the integral is 
controlled by £'k(S+_). By a Hardy estimate, we have 



The energy i?K,3 controls the integral term on the right. 

To control the second, we prove a preliminary decay result for < < t/2, 

Jo 

1/2 / rr. \ 1/2 



<ry't-l^K.3(t) + + (^K,3(0)l/' + i?T,,7(0)l 
<^-l/2+C|a| / ^^^3(0)1/2 + i?T,,7(0)l/2) . (5.2) 



/2 



In this calculation, we have used the spherical Sobolev estimate to relate the integral 
along the curve of constant (t, 9, (p) to one over the hypersurface of constant t on 
which i?K,3(i) is evaluated. This provides uniform control on the end point. Since 
we were interested in the integral over the sphere, it was not necessary to introduce 
the two, additional derivatives, so that 

aV^2|^|2^3^_^ 



52 Jr 



j^^ ^ a^lC+^l^dV* + (^K,3(0)^/^ + i?T„7(0)i/^) . 
itly small |a|, 

^,k(S+ ) <a2(^K(S+ ) + (^K,3(0)i/^ + i?T,j(0)i/^)). 



'52 Jr. 

Thus, for sufficiently small \a 



Theorem 4.13 gives us a bound on the (K,(7k) energy, 

^(K,,,)(S+ )<«^''"'^(k,,k)(So)- 

Since the (/k energy is small relative to the K energy and the initial data, we have 
the same estimate on the K energy. □ 

This control of the energy allows us to integrate along almost null surfaces to 
control the solution ip near null infinity. 

Theorem 5.4. There are constants a, C , and C , such that if \a\ < a and ip is a 
solution to the wave equation Dip = 0, then for all t > 0, r > r^, iO,<P) G "S*^, if 

u- > 0, 

and if w_ < 0, 
m,r,9,P)\ 

< Cmax{l, ~u_}-'/\-' (^K,3(0)^/' + i?T,,7(0)^/^ + supj:, (r^/^^)) . 
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Figure 2. Approximate illustration of a surface of constant u_ 
through the point {t, r*) and extending back to the surface t' = 2r^. 
Angular variables have been suppressed. 



Proof. As usual, we will start by working with the transformed function ijj which 
satisfies (4.1 ). 

With the control from the previous lemma, we can integrate along null curves. 
Consider the curve, C, from {t,r^,9,(f>) with tangent ^+ and parameter A. Along 
this curve, since the component of ^_ vanishes so rapidly, A will be uniformally 
equivalent to t, to r*, and to w+. 

First, we consider when u_ > 0. When the curve reaches the surface — t' /2, 
since w_ = t - n + 0(1), it will reach a point <, 6*', (/)') = ((2/3)(i - n) + 
0(l),(l/3)(t-n)+0(l),6',(/)+0(l)), at which will be <-n + 0(l) = u_ + 0(l). 

Integrating along the curve, 

+ o{i),^u^ + o{i),e,ci> + o{i)\ , 



1/2 



1/2 



lu^ + o{i),lu^ + o{i),e,(j) + o{i) 



^++o(i) 
-+o(i) 



1/2 



+ iii_ + o(i)ri+^i'^i^K,3(o)i/' 



1/2 



Applying second-order symmetries, integrating in the angular variables, and apply- 
ing a spherical Sobolev estimate in the usual way allows us to replace the remaining 
integral in this expression by the K energy. The end point is controlled by the 
preliminary estimate (5.2) in the previous lemma. Combining these we have 

W*, r*,^, 0)1 <(-(«+ + 0(1))-^ + (^/-+ 0(1))-^)'^' ^K.3(S+)i/2 

+ «;' + ^"" (^K,3(0)l/^ + i?T,.7(0)l/^) 

\ -1/2 

^M+(max{l, / 
To obtain the desired estimate, we note that ^ r~^|V'|- 



< max{l, u_} 



C\a\ 



/2 
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When ii_ < 0, a similar argument applies, except that there is no loss 
arising since the hypersurface of integration lies under the hypersurface u_ = and 
except that the endpoint in the integration is at (0, r^, — t+ 0{1),6, (j>), so that the 
endpoint is estimated by 

t^r., + 0(1), 6*, 0)| <(r, - t + 0(1))-^/^ sup(r-3/2^). 

This dictates the decay rate, since it is slower than the decay rate from — 
u_)/(w+(max{l, —u-})). Thus, the second estimate of this theorem holds. □ 

Appendix A. Relevance and global structure of Kerr 

In this section, we review the physical relevance and global causal structure of 
the Kerr spacetime. Most of this description can be found in most introductory 
relativity texts (i.e. [551 US])- 

The Lorentz metric for the Kerr spacetime is most simply given in terms of the 



Boyer-Lindquist coordinates {t,r,9,(f>) by (1.1 1. For a fixed constant r+, one 
might imagine the region (i, r,6,(f)) g M x (_R, oo) x to represent the region of a 
spacetime which is a vacuum outside some astrophysical object of radius less than 
or equal to The set {t, r,6,(j)) e M x (i?, oo) x S'^ can be endowed both with the 



Kerr metric in Boyer-Lindquist coordinates ( 1.1 ) and with the flat Minkowski metric 
in spherical coordinates 77 = — di^ + dr^ + r^(d9^ + sin^ 6'd0^). In this coordinate 
patch, as r — > 00, all components of g approach those of rj, which suggests that the 
Kerr metric should be seen as asymptotically flatj^ 

At this point, it is useful to review the structure of infinity for the Minkowski 
spacetime. The Minkowski spacetime can be rewritten in terms of the coordinates 
(C/+, 9, (f) where U± = arctan(t±r). If one then multiplies the Minkowski metric 
by the conformal factor cos^ Uj^ cos^ U- , one finds that the Minkowski spacetime 
can be conformally embedded into a compact subset of a larger spacetime [28 . The 
boundary of the embedded region is called the conformal boundary or points at 
infinity. Since conformal transformations preserve null geodesies, it is relatively easy 
to see that all geodesies terminate (as t 00) on the hypersurface U+ — ■n 12 and 
originate (as t —oo) on the hypersurface C/_ — —tt/2. These surfaces are referred 
to as future and past null infinity, and . All timelike geodesies terminate (as 
t 00) at future timelike infinity they originiate (as t —00) at past timelike 
infinity i~ , and all spacelike geodesies originate and terminate (both as r — > 00) at 
a single point, spacelike infinity, zq. Since the Kerr spacetime is asymptotically fiat, 
it is possible to construct the hypersurfaces X+ and X~ on which out-going (with 
t — > 00 and r ~* 00) and in-going (with t —>■ —00 and r — > 00) geodesies terminate 
and originate respectively. 

The Kerr spacetime, and especially the subcase of the Schwarzschild spacetime 
where a — 0, provides the most important illustrations of the concepts of asymptotic 
fiatness, black holes, and event horizons, which are central to general relativity. In 
an asymptotically fiat spacetime, a black hole is the complement of the past of 
future null infinity, and its boundary is the future part of the event horizon. The 
event horizon must be a null hypersurface. The domain of outer communication, or 
more simply the exterior region, is the intersection of the past of future null infinity 
and the future of past null infinity. Ignoring the physical interpretation given above, 
one can continue the Kerr spacetime as a vacuum solution of Einstein's equations by 
continuing it to r < i? and even to r < When \a\ < M , the maximally extended 



^In spherical symmetry, a = 0, it is easy to find explicit, globally smooth solutions of Einstein's 
equation coupled to matter models which can be written in spherical coordinates, in which there 
is a spherical material body from r = to r = R, which contain vacuum for r > R, and in which 
the metric exactly coincides with for r > R. We are not aware of such solutions for a =^ 0. 

^We omit the precise definitions of asymptotically flat|43| and asymptotically simple [28l. since 
they are quite technical. 
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Kerr spacetime contains a black hole, the exterior region is given by r > r_|_, and 
the event horizon is given by r = r_|_ . 

The parameters M and a can be understood by comparing the Kerr spacetime 
in Boyer-Lindquist coordinates to the Minkowski metric in spherical coordinates. 
Timelike geodesies, representing the paths of small physical bodies, in the Kerr 
spacetime do not follow the same trajectories as geodesies in the Minkowski metric. 
The deviation between these geodesies can be treated as a relative acceleration, and, 
as r ^ oo, this acceleration corresponds to that generated by the gravitational force 
from Newtonian mechanics (also in spherical coordinates) generated by a central 
object with mass M. This is one reason why the parameter M is interpreted as the 
mass of the black hole, even though the Kerr spacetime is a solution of Einstein's 
equation in vacuum. 

Alternatively, Einstein's equation can be treated as a Lagrangian theory [SJ, and 
the parameters M and o can be understood as energies associated with a vector 
fields. Normally, an energy is computed on a hypersurface, but, under certain 
conditions, one can apply integration by parts so that the energy can be computed 
as an integral over the boundary of the hypersurface. For any asymptotically flat 
spacetime, this allows energies to be defined with respect to vector fields defined 
near infinity in Minkowski space. In Minkowski space, the energy associated with dt 
is referred to as simply the energy. In relativity, one expects that, with the correct 
choice of units, the energy and mass coincide. In Minkowski space, the energy 
associated with a rotation about an axis is called the angular momentum about 
that axis. Computing these quantities on two-surfaces very close to zq, one finds 
that the mass of the Kerr black hole is Af , that the angular momentum about the 
axis of symmetry is Ma, and that the angular momentum about orthogonal axes is 
zero. 



The metric (1.1) is clearly smooth in the exterior region, since neither A nor E 



vanish. The roots of A are at 

r± = M ± -a2, 

and the roots of S occur when both r and a cos 9 vanish. Thus, when < |a| < 
M, the metric is clearly smooth where r > r+, where r_ < r < r+, and where 
r < r_ except at (r, 0) = (0,7r/2). We shall refer to these regions as the exterior, 
intermediate, and deep-interior regions respectively and as block I, block II, and 
block III respectively. Using the coordinates (u_, r, 0, 0,^) (with u_ defined in the 
introduction, r and the standard Boyer-Lindquist coordinates, and 0* = — < 



given in section 5.2 1, it is possible to construct a smooth coordinate chart which 
covers a copy of both block I and block II. 

Through similar choices of changes of coordinates, various copies of the blocks 
can be be smoothly joined together to form the maximal extension of the Kerr 
manifold as illustrated in figure |3] The copies of block II and block III meet where 
r = r_ . The two-dimensional surfaces where two copies of block I and two copies 
of block II meet is also a regular point of the maximal extension and is called the 
bifurcation surface. The diagram in figure [3] is also a conformal diagram, which 
includes the conformal boundaries of the exterior regions and also of the deep- 
interior regions. The deep-interior regions contain singularities, in the sense that 
there are affinely parameterised geodesies along which (r, 0) — s- (0,7r/2), on which 
the affine parameter remains bounded, but which cannot be extended, and it is not 
possible to smoothly extend the spacetime so that the geodesies can be extended 
to (r, 0) — (0, 7r/2). The diagram does not reveal other patholgical behaviour in the 
deep interior, in particular, the presence of closed timelike geodesies. 

The subcritical (|a| < M) Kerr spacetime is expected to have a crucial physi- 
cal role. From the singularity theorems, it is known that under a wide range of 
circumstances, spacetimes will develop singularities. The weak cosmic censorship 
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Figure 3. The conformal diagram for the Kerr spacetime: Angu- 
lar variables have been suppressed. The boundaries of blocks are 
shown in solid lines. Boundaries at infinity and singularities are 
indicated by dotted lines. The event horizon at r r+ and X+ 
has been indicated for one exterior region. The maximal extension 
continues with this pattern infinitely forward and backward in time. 

conjecture asserts that, in asymptotically flat spacetimes, future null infinity is sep- 
arated from any (future) singularity by an event horizon, as long as some genericity 
condition holds. The statement of the genericity condition remains to be deter- 
mined. The no hair theorem, or Kerr uniqueness theorem, states that the Kerr 
family of spacetimes (including Minkowski as M — 0) is the unique class of station- 
ary spacetimes satisfying Einstein's equation in vacuum (see |38j for early work, 
and for recent progress) . It is common in physics that systems will relax towards 
stationary configurations, and, in particular, it is widely expected that solutions 
to Einstein's equation containing a single black hole will asymptotically approach 
a Kerr solution plus gravitational radiation going to X"*", at least in the exterior 
region. A preliminary step in showing that this is true would be to show that the 
Kerr spacetimes are dynamically stable, i.e. that a small perturbation of the initial 
data which leads to a Kerr spacetime will generate a solution to Einstein's equation 
which asymptotically approaches a Kerr spacetime plus radiation. 
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